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ABSTRACT 

The polarization transfer coefficients of a relativistic magnetized plasma are derived. 
These results apply to any momentum distribution function of the particles, isotropic 
or anisotropic. Particles interact with the radiation either in a non-resonant mode 
when the frequency of the radiation exceeds their characteristic synchrotron emission 
frequency, or quasi-resonantly otherwise. These two classes of particles contribute dif- 
ferently to the polarization transfer coefficients. For a given frequency, this dichotomy 
corresponds to a regime change in the dependence of the transfer coefficients on the pa- 
rameters of the particle's population, since these parameters control the relative weight 
of the contribution of each class of particles. Our results apply to either regimes as well 
as the intermediate one. The derivation of the transfer coefficients involves an exact 
expression of the conductivity tensor of the relativistic magnetized plasma that has 
not been used hitherto in this context. Suitable expansions valid at frequencies much 
larger than the cyclotron frequency allow us to analytically perform the summation 
over all resonances at high harmonics of the relativistic gyrofrequency. 

The transfer coefficients are represented in the form of two- variable integrals that 
can be conveniently computed for any set of parameters by using Olver's expansion of 
high-order Bessel functions. We particularize our results to a number of distribution 
functions, isotropic, thermal or power-law, with different multipolar anisotropies of low 
order, or strongly beamed. Specifically, earlier exact results for thermal distributions 
are recovered. For isotropic distributions, the Faraday coefficients are expressed in the 
form of a one-variable quadrature over energy, for which we provide the kernels in 
the high-frequency limit and in the asymptotic low-frequency limit. An interpolation 
formula extending over the full energy range is proposed for these kernels. A similar 
reduction to a one-variable quadrature over energy is derived at high-frequency for a 
large class of anisotropic distribution functions that may form a basis on which any 
smoothly anisotropic distribution could be expanded. 
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1 INTRODUCTION AND MOTIVATION 



The transfer of the intensity and polarization of radiation propagating in a relativistic magnetized plasma has been a 
long standing problem which still remains i ncompletely solved. Its potential diagnostic applic ations have long been known 



JCioffl fc Jone: 



Thiebaut et al 



1980] ; I Jones fc O'Delll 119771 ) and it has recently attracted renewed interest (jBrentiens fc de Bruvnl 12005 



bold ) in connection with the study of the magnetic field structure of synchrotron-emitting objects where 



the high-energy emitting particles may constitute the bulk of the plasma. Polarization data constrain the structure of geo- 
metrically complex sourc es and of their magnetospheres. Astrophysical environments such as t he vicinity of the ce ntral black 



hole in t he Milky Way (lAgolll200d : Bcckcrt 



2003 



Shchcrbakov 



2011 



jet-sources, galactic Jstirling et aL 2004 ]_or extra 



galactic Jfioman "fc Wardkfaod ). pulsars jKanbach et allbooj Dyks et allbood : IPetri fc Kirkll2005l : iMc Donald et alj|20icl : 



Yuen et al]|2012l ) or disks around black holes jpovciak et aljboosl ) havebeen studied by these methods. The perspective s 



of using polarization data to constrain magnetic fields have improved strongly by the deployment of LOFAR ( Beck 2009) 
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and will blossom with the advent of SK A ( Dewdnev et al. 20091; Carillil 2004 ) , ranging from the largest scales ( Gaensler et al 
20041 ; iFeretti fc Johnston- Hollittll20o4 ) to the smallest ones (|Bicknelj|20o4 ). 



Undestanding polarization data rests on the availability of good polarization tranfer codes and on the knowledge of the 
transfer coefficients that would modify the intensity and the polarization of light as it propagates in the medium. Synchrotron 
radiation plays an important role in the emission of magnetized high energy sources. The polarization properties of its emission 
and absorption caused by the high-energy population are well known. The non-dissipative transfer coefficients of Faraday 
rotation and Faraday conversion from circular to linear polarization have been more difficult to derive. Substancial progress 
has nevertheless been achieved in time (|Trubmkovlll958l : IZheIezmakov|[l9r^ 

Huang fc Shcherbakovll201lh . In this paper we buildup on this research by deriving explicit expressions for the polarization 



transfer coefficients in a relativistic homogeneous magnetized plasma for any isotropic distribution function and for a large 
class of non-isotropic ones. These results are obtained by a new exact method of derivation of the high frequency conductivity 
tensor of a magnetized relativistic plasma which involves, at some late stage of the derivation, well-justified approximations. 
Our results are shown to exactly coincide with known exact results in the case of thermal distributions. 

Section [2] is devoted to introducing the polarization transfer equations of uniformly magnetized medium and deriving ab 
initio the corresponding conductivity. Section presents in anticipation the end-results of our derivation of the polarization 
transfer coefficients as a double regular quadrature and illustrates it for isotropic and anisotropic distributions, where it can 
sometimes be reduced to a single quadrature. This unconventional presentation has been chosen because the derivation of 
these results from those in Section [2l which is carried in Section [4] and Section [5l is technical in some respects. These two 
sections can be skipped accordingly on a first read. Section [6] focuses on analytical results applicable to isotropic distribution 
functions in the low-frequency limit. Section wraps up. 



2 FORMULATION 



2.1 Description of polarization transfer 

The polarization properties of radiation are represented by a column vector of four Stokes parameters, I = T (I, Q, U, V), T M 
denoting the transpose of the matrix M. The Stokes vector I depends on the frequency and direction of propagation of the 
radiation, as well as on the point in space and on the time at which it is observed. In the limit where the frequency is higher 
than any gyrofrequency, the Stokes vector satisfies a radiative transfer equation in which the generalized emission coefficient 
is a four-component column vector W = t (Wi,Wq,Wu,Wv) and the generalized absorption coefficient is a 4x4 transfer 
matrix. This transfer equation can be written, for radiation of frequency uj and group velocity (based on an isotropic mean 
conductivity) v s as: 
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The four matrix elements Kn, Kiq, Km, Kiv an d the emission column vector W result from dissipative effects. The emission 
column vector W represents the emission rate and polarization of the spontaneously emitted radiation. The elements Kqu, 
Kqv, Kuv of the transfer matrix result from non-dissipative effects caused by the slowly varying phase difference between 
different components of the electric vector of the radiation. The elements of the matrix in equation {TJ can be expressed in terms 
of the components of the h igh-frequency conductivity tensor er(k, ui). Equation CD has been establi shed by a number of authors 



and by different methods ( Sazonov fc Tsvtovitch 19681 ; Hevvaertsl 19691 ; Zheleznvakov et al. 19741 ). A direct derivation of the 



absor ption-like term is presented below. The derivation of the emission term is less straightforward (|Sazonov 
196Sl ). 



Tsvtovitch 



Since the wavelengths of the plasma eigenmodes of a given frequency u do not differ by much from each other at high 
frequency, their phase difference evolves slowly in space and time. As a result, the correlations between these phases, which 
determine the state of polarization, are maintained over a length much longer than the wavelength, justifying a radiative 
transfer approach. When the plasma frequency lo p exceeds by much the gyrofrequencies Qa of all particle species A, the 
wavelengths of the eigenmodes still differ by a little amount when ui 3> although the dispersion relation may depart 
from the vacuum one when ui is of order w p . In any other case, the difference of the phases of the transverse eigenmodes 
evolves on a wavelength scale and the radiative transfer approach is inappropriate. Assuming there are only two transverse 
eigenmodes, if their phases would be correlated at some point and time, the rapid growth of their difference would reduce 
this correlation to nothing in only a few wavelengths. Thus, the polarization state of radiation supported by modes of largely 
different wavelengths for a given frequency reduces to that of the sum of two mutually incoherent eigenmodes. Equation CD 
is then only relevant in the high frequency limit, u ^> Qa- 
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We now show how the elements of the absorption matrix relate to those of the conductivity tensor. We start from the 
propagation equation of the electric field E(r, t), deduced from Maxwell's equations: 

iw- aE + -"f -"• < 2 > 

Our equations are written in the MKSA system of units, where eo is the dielectric permittivity of vacuum and fio its magnetic 
permeability. The velocity of light in vacuo is c = (eo^o)~ 1 ^ 2 ■ In equation p)l. E is the electric field and J the current density. 
The electric field is assumed to be transverse, so that divE = 0, which is a good approximation at high frequencies. The 
microscopic electric field is then Fourier expanded, its complex vectorial amplitude Ek, u being regarded as slowly depending 
on space and time, which is equivalent to coarse-graining in Fourier space. The electric current density J(r, t) is similarly 
expanded. The phase, and possibly the amplitude, of the microscopic electric field is a random variable. Omitting the indices 
k, oj for simplicity, each term of the Fourier expansion can be written as 

E(r,t) = E(r,f) e l(k,r ^ 4) . (3) 

The current in equation ((2} is the sum of an induced current J', which results from the electric field E(r, t) and a spontaneous 
current J , the microcurrent created by the individual plasma particles, which gives rise to spontaneous emission, while the 
induced current gives rise to the absorption- like term in equation ^ . This latter part of the current is given in terms of the 
electric field by a linear non-local and causal conductivity tensor operator: 

J'(r,t) = j^ dTjjjd 3 R S(R,-r)-E(r-R,t-r). (4) 

The conductivity operator can be split into a scalar isotropic part, S s , and a tensorial anisotropic part Si. The isotropic 
part may, for example, represent the dispersive properties of a cold population of relatively large density, neglecting the 
magnetization of these particles. An isotropic conductivity has no effect on the polarization but affects the propagation 
properties of electromagnetic waves. No component of the conductivity need to be separated out when the propagation is 
represented well enough by the vacuum dispersion relation. The tensorial anisotropic part consists of the anisotropic residual 
of the conductivity, not included in S s , and of the conductivity of populations which negligibly contribute to the dispersive 
properties of the plasma waves. The frequency u) supposedly being much larger than all gyrofrequencies, the component Si 
must be small. Its contribution to the dispersion properties may be neglected but it nevertheless contributes to the polarization 
transfer matrix. 

The form of the electric field in equation Q is inserted in equations ([2]) and Q. When acting on the imaginary exponential 
phase factor the operators dt and V are of order uj and k respectively, but when acting on the slowly varying amplitude E(r, t) 
they are of order 1/T and 1/L, T and L being the time and length scales on which E varies, which are much longer than the 
wave period and wavelength. Only the terms of first order in l/(ioT) and \/(kL) are retained. For consistency, this requires 
that E(r — R, t — r) in equation Q be expanded as: 

E(r - R, t - r) = E(r, t) - r d t E - (R ■ V) E . (5) 

The time and space integrations in equation @ then involve the Fourier transforms of S B and Si, cr s (k, oj) and <r(k, ui) (no 
index 1), the latter being a tensor, and the derivatives of <r s with respect to lo and k. This results in: 

^ _ C V + — a<\ E + (2110 — l^A) 9t E + ( ( 2zc 2 k + -V k a s ) • v) E = -™ <r(k, u ) ■ E . (6) 

eo J \ eo du J \\ eo J J eo 

The term of dominant order in equation © is the first one. The two other terms on the left are of first order in l/(uT) and 
l/(kL) and the term on the right is also regarded as small. Equation © is satisfied at the dominant order when the factor 
of E vanishes, which means that the electric field consists of fluctations, the frequency of which is uik, the solution of the 
isotropic dispersion relation: 

uj 2 — c 2 k 2 H os(k, ui) = . (7) 

eo 

The corresponding group velocity, v 9 , is obtained by differentiating equation (0 with respect to uj and k. When a s represents 
the dispersion due to a cold unmagnetized plasma, which we assume for simplicity, uja a does not depend on to nor on k. With 
these simplifications, the three sub-dominant terms in equation ((6]) reduce to: 

(0 t + (v. • V)) E = - - ^ , (8) 

where the amplitude E, which slowly varies with position and time, is only non-vanishing when the wave vector k and the 
frequency uik are linked by the dispersion relation in equation (0. An equation for the correlation tensor of the electric field 
of this radiation, (E ® E*} can be derived from equation ((Sj by tensorial multiplication by E*, the superscript * denoting 
complex conjugation, and averaging over the statistical distribution of the phases (and possibly of the amplitudes) of the 
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components of E. This operation is represented by brackets. The electric field supposedly being transverse, it pertains the 
plane perpendicular to the propagation direction n and may be represented by its components on two abitrarily chosen unit 
basis vectors ei and e2 in this plane. An intensity tensor with components iy relative to this basis that is proportional to 
the two-dimensional electric field correlation tensor must be defined such that it has the dimension of a specific intensity. The 
Stokes parameters of the radiation are related to this intensity tensor by 



k 2 dk 
du) k 



eo{E 2 El) 



so(EiE%) 
soiEzEZ) 



I + Q U + iV 
U-iV I-Q 



(0) 



The definition of the Stokes paramet er V is a matter of convention which fixes its sign, since (-Ei-Ef) could have been taken 
to be proportional to U — iV instead. Shcherbakov fc Huan j ( 2011 ') use a convention opposite to ours. The transfer equation 
for the components ly can be deduced from equation J8} which writes, using the dummy index rule 



(ft + (v a • V)) U 



— {oipSjq + Si P a* q ) 



(10) 



By equation ([9|, the term on the right of equation (|10|l can be be converted into the transfer matrix term in equation (fl}, some 
elements of which turn out to involve the hermitian (dissipative) part of the conductivity while others involve its anthermitian 
(non-dissipative) part. These parts are respectively defined by: 



h 1 / . * \ a 1 / * \ 

a H = 2 \ aii + a i i > ' ff « = 2 ^ j ~ ji > ' 

The elements of the transfer matrix in equation fT]) are given in terms of these parts of the conductivity by: 
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Eo ^£0 £o 

These matrix elements refer to the components of the local conductivity for the wave vector k and frequency cjk related to it by 
the dispersion relation in equation Q. The relations in equation (1 1211 between the transfer coefficients and the elements of the 
conductivity matrix agree with those in Ishcherbakov fc Huang ( 2011 ) given the relation between conductivity and dielectric 
tensors for perturbations varying as in equation Q and accounting for the fact, mentioned above, that their definition of the 

Stokes parameter V differs from ours by a sign. 

The polarization properties of the synchrotron e mission are well known, both in vacuo ( Westfold 19591 ) and in a cold 
plasma 1 Tsvtovitch 1950l : Razin 1960l ; Ramatv 19681 ). The dissipative absorption coefficients Ku, Kiq, Kiu, Kiv given 
in equatio n (I12H in terms of the components of the he rmitian part of the conductivity have also been presented in the 
literature ( Ginzburg fc Sirovatskii 19691 : Sazonov 1969a ). These components can be calculated using the same well-known 
approximations which also provide the emission co efficient. The dis sipative absorption can be d erived from the emission by 
Einstein's coefficients method, ignoring dispersion (jwild et alJll963h , or not (|ZhelezniakovHl967l ). 

The non-dissipative coefficients Kqu, Kuv, Kqv have proved more difficult to calculate. This paper concentrates on 
their calculation in the case of ultrarelativistic plasma particles immersed in a static uniform magnetic field Bo. The motion 
of the particles is described in the Vlasov approximation. We neglect the dispersive properties of a cold population that might 
be present, so that the dispersion relation (J7) supposedly reduces to uj = ck. 



2.2 Formal solution for the high frequency conductivity tensor 

The conductivity tensor is found by calculating the high frequency current J resulting from a high frequency electric field E 
by solving the linearized Vlasov equation for the perturbation /(r, p) of the distribution function, 

dt f + v ■ V/ + (qv x Bo) ■ V p / = -q{E + v x B) • V p / , (13) 

where q is the charge (with its sign) of the particle species considered. In equation (|13|) fo (p) is the unperturbed homogeneous 
distribution function of the species considered, that only depends on the component p« of the particle's momentum along the 
magnetic field and on the modulus p± of its component perpendicular to it, normalized to the density of particles, so that 

2TTp ± dp ± dp n f (p ± ,p n ) =n, (14) 

where n is the volume density of this species of particles. The gradient with respect to position is denoted as V and the 
gradient with respect to some other vectorial variable u, such as the momentum p, is denoted by V u - The first three terms 
of equation (|13[) represent the time-derivative of / following the unperturbed particle's motion, described by the position and 
momentum r(t') and p(t') of this particle at time t' , following the unperturbed motion. At time t, the particle is at r with 
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momentum p. We then obtain the perturbation / by integrating in time the right hand side term of equation (|13|) following 
the unperturbed particle's motion, which gives 



f(r,p,t) = -q f «ft'(E(r(0,0+v(t')xB(r(0,*')-Vp/o(p(t'))- 

J — oo 



(15) 



The integration is over all times t' earlier than t. The unperturbed motion is easily expressed in a frame where Bo is along 
the Z axis so that Bo = B^ez and the two unit vectors ex and ey are perpendicular to it, but otherwise unspecified. The 
velocity of a particle can be written as 

v = v± cos 4> ex + v± sin (f> ey + Vu ez , (16) 

where v± and vu are conserved by the unperturbed motion, as are also the modulus p of the momentum, the associated Lorentz 
factor 7 and the particle's pitch angle #. The angle <j), the gyration angle of the particle, rotates in time at the synchrotron 
frequency Q t , which for non-relativistic particles reduces to the cyclotron frequency Q. These frequencies, which have the sign 
opposite to that of the charge, s q , depend on the rest mass m of the particles and are given by 

f2* = — - ^° . Q — — - — - , and s„ = sign(g) . (17) 

7m m 

At a time t' = t — r earlier than t, the components of the velocity of a freely-moving particle having at time t a gyration angle 
<j> and a position X, Y , Z were 

vx(t — t) = v±(t) cos(0 — , Vy{t — t) = v±(t) sm(<fr — J1*t) , vz(t — t) = v^(t) . (18) 

The position of that particle at time t — r was 

X(t-T) = X-^ (sin0-sin(0-ft*-r)), Y(t - r) = Y + ^ (cos - cos(<?!> - 0,t)) , Z(t - r) = Z - V\\ r. (19) 

The magnetic perturbation is related to the electric one by the Faraday equation, curlE = — dtH. The perturbation in the 
distribution function which develops as a result of an electric field perturbation such as that in equation ([3j can be written 
as /(r, p, t) — /(p) exp(i(k • r — tot). It depends linearly on the perturbation field by equation (| 15p and can be expressed in 
terms of operators acting on the unperturbed distribution function, such as the anisotropy operator D, defined by: 

D(fo) vanishes for an isotropic distribution function. The wave vector k is taken to be in the X — Z plane and has components 
kx = fcsina and kz — fccosa, the propagation angle a being comprised between and n. When fully expanded, the solution 
for / can be written from equation (|15|l as 

/(p) = — — / dr exp ( i ( (u — k\\V\\ ) r — —(sin <fi — sin(</> — Q*t) 
U Jn V V S2* 



+ fc||D(/o)^ + E z (oj^ - cos(9i - fi*T) fc_L-D(/ ) 



(21) 



The associated Fourier coefficient of the current density is the average value of gv wheighted by /(p). Using the solution for 
/ in equation (JUJ), each component of the current appears to be a linear function of the electric field components Ex, Ey 
and Ez from which the components of the conductivity tensor may be identified. Each one involves a fourfold integral on 
the variables p«, p±, <f> of the particle momentum and on the delay time r. As can be seen in equation \21\ , integration over 
the delay time involves imaginary exponentials with an argument linear and trigonometric in r. We refer to these integrals 
as phase integrals. Expanding the trigonometric functions of (</> — Q»t) on the second line of equation (|21|) in imaginary 
exponentials, three phase integrals appear, that can be written as 



Pe(u), k±, fc||, v) = J dr exp \ i\(<*J — fc||t>||)r — (sin^ ~ sin(<p — O^t)) + e{cj> — 0„,t)^ 



(22) 



The index e may take the three values (—1, 0, +1), or (— , 0, +). At this point, it is necessary to take care of the sign of the 
charge, and to introduce dimensionless variables, simpler than those in equation (|22p . by: 

kj_v± w-fc||i>|| u) ck\\-v\\k 

X = W a = ^hT' U = W\' ff0 = tt8mQ ' m = (23) 

The variables x (not a coordinate), a and u are all positive and usually large since u) = ck 3>| |. The variable zu, which 
plays the role of an angular variable, may be positive or negative. The modulus of w assumes values of the same order as a. 
The properties of the variables vo and a, which may replace p\\ and p± as dynamical variables of freely moving particles, are 
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detailed in Appendix [X] From equation (|21[) . the current density can be written in terms of the three phase integrals as: 
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(24) 



3 RESULTS 

We derive in section [4] the components of the conductivity tensor for a relativistic plasma at frequencies much higher than 
the gyrofrequency and carry out two out of the four integrals involved in equations (|22|) and (|24|) . In particular cases, the 
corresponding quadrature may be pursued one step farther. 

In this section we anticipate the polarization transfer coefficients and illustrate them for isotropic and anisotropic distri- 
butions. 



3.1 Dissipation-less polarization transfer coefficients 



3.1.1 Quasi exact expressions 

We obtain in section l4~4l a general expression for the Faraday rotation and conversion coefficients / and h, defined / = Kqu/c 
and h = Kuv/c, where Kqu and Kuv are given by equation (|12|) . These coefficients can be explicitely written as in equations 
(18611 (I87p . which we reproduce here neglecting the principal value terms that are shown in section |4~41 to be safely negligible. 
These equations are therefore almost exact. They can be written as 



f 



,2 



m 3 c 3 , , dFo 
— 5— dzu da tux — 
sin a da 



j' a {x)N a {x) + 
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2 J a {x)N a {x)) + 
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(25) 
(26) 



The particle's dynamical variables zu and a, defined in equation (|23p . incorporate properties of the radiation of interest. Much 
of the dependence of / and h on the propagation angle a and on the frequency is in fact hidden in these variables. The 
distribution function F is normalized to unity (in the sense of equation (I75[) below), the particle density n being implicit 



in the square of the specie's plasma frequency cu pl — nq /meo- The integrals in equations (|25[) and (|26[) involve two- variable 
kernels in which one of the variables, a, appears as a continuous Bessel function index. This index is larger than their argument 
x, as can be seen from equation (|23|l . but possibly not much larger, which complicates the search for suitable approximations. 

There are essentially two very different regimes, defined in section [S] for the wave-particle interaction which we refer to as 
non-resonant (NR) and quasi-resonant (QR). Different particles make very different contributions to the transfer coefficients 
depending on whether they interact in the NR or in the QR regime. Equations (|25|) and (|26[) however encompass all regimes, 
being written in a form that does not require any regioning of the zu-a domain for their evaluation. The domain of validity 
of specific approximations, non-resonant or quasi-resonant, to the functions in equations p5[) and (|26[) are different and 
complementary (figure |SJ) . 



3.1.2 High-frequency limit (HF) 

It is shown in section [BTTI that when the inequality 

372sin " 1 " 1 < 1 (27) 

id 

holds true, only the non-resonant domain, defined in that section, contributes to the transfer coefficients. For a given value of 
the Lorentz factor 7, the inequality (|27[) is satisfied when the radiation's frequency ui exceeds the characteristic synchrotron 
frequency w c (7, a) for particles of this energy travelling in the direction of the radiation (equation (flOT}). We therefore refer 
to this situation as the high-frequency limit. When all particles in the distribution interact with the wave in this regime, the 
integration on the non-resonant domain actually extends over the full zu-a domain. In this case, the Faraday coefficients / 
and h result from equations (|115[) and (|119[) . expanded to the first non- vanishing order in Q/w: 

^ da ( 2 (o^-xiy/i ) (28) 
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Equations (|28[) and (|29[) do not assume isotropy of the distribution function Fo but are only valid when the quasi-resonant 
contribution can be neglected, which makes them less general than equations ([25} -([26]). 



3.1.3 Medium- frequency (MF): cutting through the quasi-resonant domain 



The integration over the non-resonant domain cannot be extended to the full domain when the inequality in equation (|27|) is 
not satisfied for a non-negligible number of particles in the distribution. This arises in particular when the parameter on the left 
of equation (|27|) is of order unity, which we refer to as a medium-frequency situation. A slightly more sophisticated approach 
is then needed, which we now describe. The double quadratures in equations (|28l) - (|29[) should in this case be extended over 
the non-resonant (NR) domain only and the contribution of the quasi-resonant (QR) domain should be added to it. When 
it is not dominant, the integral over the QR domain may be approximately evaluated. Our proposed approximation rests on 
separating the non-resonant from the quasi-resonant domain, and expressing the integral over the quasi-resonant one by using 
a very simple approximation. The boundary Sqr between the NR and QR regions is where the variable g defined in equation 
(|94[) below is near unity. Since x ~ a in the vicinity of this limit, Z3qr is represented in the zu-a plane by the equation 



Bqr(zU,o) = ZU 



d2/3 4/3 



+ a 2 = . 



(30) 



The sign of Bqu(zu ,a) tells whether the point zu-a is in the non-resonant domain (-Bqr > 0) or in the quasi-resonant one 
(Bqr < 0). The QR domain is then cut out of the zu-a domain, the remainder constituting the properly-defined NR domain. 
The kernels of the integrands in equations (|28[) - (|29|) are kept at their non-resonant values when in the NR domain. For perfect 
accuracy, the kernels in the QR region should be taken to be those in equations (|99p - (|lU0|) . When however the QR region does 
not contribute predominantly, these complicated kernels can be replaced by simple approximations. Their effect essentially 
being to limit the growth of the modulus of the non-resonant kernels when the QR region is reached, we replace them by linear 
interpolations in zu, at fixed a, between the values of the non-resonant kernels evaluated on the border Sqr for that value of 
a. Let On(y) denote the Heaviside function, equal to unity for positive argument y and zero otherwise. The above-described 
approximations to / and h can be written as: 



/mf = 2ns q r- 



- dzu da zu 



(-) 



+ e H 



■ dzu da 



H 2 (a 2 - x 2 ) 3 / 2 ' ~ a 2 (a 2 
\a 2 -x 2 ) + a 2 ) x 2 (4a 2 +x 2 ) 



T 2 

^QR 



8 (a 2 - x 2 ) 7 / 2 



4r) 3/ 
e H 



2zqr 



dFp 
da 

i° 2 - 



(31) 



2 \ , 2 2 
' ZQRj + a X QR 



(4<7 2 +z|r) 



8 (a 2 



QR 



y/2 



dF 
da 



(32) 



where <d^ = 0h(±-Bqr(^, a)) and zuQn(a) and iqr(<t) are the values of \zu\ and x on the boundary Sqr for the given value 
of a (equation (|30p ). Equations (|31[) and (|32[) provide approximate values of the transfer coefficients / and h valid for isotropic 
as well as for non-isotropic distribution functions. The approximations in equations (|31[) and (|32|) will be illustrated on the 
example of a thermal distribution function in section 13.2.41 



3.1.4 Low-frequency limit(LF) 

When the reverse inequality to equation (|27[) applies in a strong sense, both the non-resonant and quasi-resonant domains 
contribute to the transfer coefficients. For a given frequency and a given propagation angle, a, this essentially happens when 
the dynamical variable a of the particle is above the limit defined by equation (|30[) for zu — or when its Lorentz factor 
7 largely exceeds the threshold defined by the reverse of inequality (|27[) . which happens when the frequency ui is much 
less than the critical synchrotron frequency, for a particle of Lorentz factor 7 travelling in the direction of the considered 
radiation (equation fll01| |). We refer to this situation as the low- frequency limit. The low- frequency contribution to the Faraday 
coefficients comprises both NR and QR parts, which we calculte for isotropic distributions in sections 16.31 and 16.41 



3.2 Kernels of polarization transfer coefficients for isotropic distribution functions 

Let us illustrate how to use equations ([25} and (|26[) on isotropic models of the distribution function. When the distribution 
function is isotropic, it depends only on the Lorentz factor 7 of the particles, but nevertheless depends on both variables 
zu and a. When dealing with isotropic distribution functions, it is therefore wise to change from the variables zu, a to the 
variables zu, 7. Integrating over zu, which then does not involve the distribution function, equations (|25[) and (|26[) can be 
written in the general form 

. dFo „iso, \ , ^p r fi 2 3 3 f°° dFo „iso/ x , 00 \ 

d<y — F (7), h^—^mcj^ dl — H (7), (33) 

where F 1 s °(-y) and ff lso (7) are one- variable kernels for isotropic distribution functions. 



U) pr \l 33/ 

/ = 3- m c 
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3.2.1 Isotropic distributions in the high-frequency limit 

In the high-frequency limit, the integration over zu in equations (|28[) — <|29f) . which defines these kernels, extends, for a given 
value of 7, over the full domain. The kernels _F lso and H 1B ° can then be written in this case as 



2 (cr 2 - a: 2 ) 3 / 2 



; tt ( 2x A {a 2 -x 2 )+a 2 x 2 {Aa 2 +x 2 ) \ 

sm a J u ,_ \ 8 (cH — x*)'!* J 

The variables a, x and (a 2 — x 2 ) can be expressed from equation (|A2[l in terms of 7, vj and parameters in equation (|23|) as 

a = 7 u sin 2 a + w cos a , x 2 = (u 2 (7 2 sin 2 a — 1) + 271* vj cos a — uj 2 ) sin 2 a , a 2 — x 2 = uj 2 + u 2 sin 2 a . (36) 

The boundaries zo-(pi) and zu+(j) of the integration over zu at a given 7 are the abcissas of the intersection of the line of 
constant 7 with the boundary of the physical domain in the vu-a plane (figure lATjl . Since this boundary is where x = 0, 
equation (|36[l implies that 



zu± (7) = u ^7 cos a ± -\/7 2 — 1 j , (37) 

The integration over in equations (|34|) and ()35|) plays the role of an angular integration. It can be performed analytically, 
the results, valid at high frequency, eventually being 

J F , H 3 F(7)=4 7r S9 ^^^(7£( 7 )-y 7 ^T) , where £(7) = cosh"^) = In (V^l + 7) , (38) 
^ s f(7) = - ^ (7 n/t^I (27 2 - 3) + £(7)) ■ (39) 



3.^.2 Isotropic distributions in the low-frequency limit 



It comes out of the discussion in section 15.41 that the coefficient of Faraday rotation / does not depart much from its non- 
resonant approximation, although corrections are necessary for a quantitative agreement in low-frequency situations. We find 
in section [6] that the isotropic kernel for / can be written in the low- frequency limit as 

eM=-T^ k £)-H' ,40) 

By contrast, the Faraday conversion coefficient h drastically differs from its high-frequency expression when there is a quasi- 
resonant contribution. The kernel of this coefficient is derived in section 16.31 in the low-frequency limit and for an isotropic 
distribution function. It can then be written as: 

„ Sb) ^(^r( 4 M 7 « (41) 



3.2.3 Exact isotropic kernels and composite approximations to them 

The results in sections 13.2.11 and 13.2.21 apply in different regions of the energy space. For distributions with only very few 
particles in the low- frequency regime, the high-frequency approximation is sufficient. When the distribution is more extended 
in energy than the threshold 7 qr defined in equation (|98[) . the kernels should account for the QR contribution in the LF 
regime. They may then be calculated by integrating over w at fixed 7 the right hand sides of equations (|25[) - (|26[) . This 
however requires a double numerical integration. Section \6 . 5 1 proposes an approximate formula interp olating between the HF 



expres sion of these kernels, valid when 7 < 7 qr , and their asymptotic LF expressions, valid for 7 3> "iqr. lHuang fc Shcherbakov 



( 201 If ) adopt a similar approach by providing fits to their numerically computed kernels and claim similar accuracy. 

For extended energy distributions, the region at 7 ~ 7 qr should not dominantly contribute to integrals such as those in 
equation (|33[). These interpolation formulae then provide about 10% accuracy on h for power law distributions and slightly 
better accuracy on /. When on the contrary the energy range close to y qr is determinant, the medium frequency approximations 
in equations (|31[) - (|32[) could be used instead or, if accurate results are needed, the transfer coefficients would be numerically 
calculated from equations (|25[) and (|26[) by double-integrating over the variables vo and a. 
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Figure 1. Left: ratio to the high frequency expression fjjp in equation (1431 ) of resp. the numerically-computed value of the thermal 
Faraday rotation c oefficient f th for the isothermal model of equation 1142 I I. based on equation H25I) . the fit formula provided for it by 
IShcherbakml j200Sf) and the region-cut approxim ation in equation d 3 1 D . as labeled. These ratios are displayed as a function of the scaling 
parameter X introduced by IShc hcrbakov (200J, defined in equation I I47II . Right: the high frequency limit of the Faraday conversion 
coefficient h th from equation 1 143 II . the region-cut formula from equation (132 P and the almost exact value numerically computed from 
equation A25I I, as a function of the temperature, T (as labeled). In both cases, the region-cut procedure around the resonant particles 
seems to improve significantly the domain of validity of the match to the exact transfer coefficients. 



3.2.4 Illustrations and comparison to known results 



Previous approaches to polarization transfer coefficients 

The polarization transfer coefficients for a thermal distribution have been previously de rived exactly, by a method entirely 
different from ours, in a form that eventually requires only one numerical quadrature ( Trubnikov 19581 ). It is therefore of 
interest to test d ifferent approxima te results against this exact result, or against results proposed for the thermal case by 
different authors. Trubnikov ( 1958T ) considers a Jiittner distribution, as in equation (|42[) below. The conductivity is derived 
by integrating the changes of the perturbed distribution function over the delay time along unperturbed trajectories, just 
as we did in section O then taking suitable moments of the components of momentum. The exponential dependence of 
the Jiittner distribution on energy can be associated to the complex exponential that arises from the evolution over time 
of the phase functions, as in equation (|57l) . The integration over the particle's energy can then be carried out, yielding 
a result that involves a modified Bessel fun ction of the second kind of a complex argument depending on the delay time 
( Melrose! [l997al ; Ishcherbakov fc Huanell201ll). The integration over the delay time then remains to be carried out. This can 
be done by different numerical strategies (|Shcherbakov fc Huanell201ll ; iHuang fc Shcherbakov 2011), yielding exact r esults. 



Else, it can be carried out analytically in an approximate way in different lim i ts ([Melrose 1997cJ ; Shcherbakov 20081 ) or by 



using the method of stationary phase ( Melrose! 1997dl ). Huang fc Shcherbakov ( 20 111 ) have extended this approach to other 
distribution functions, the integration with respect to delay time being performed first for a given value of the Lorentz factor. 
These authors numerically calculate transfer coefficients for Dirac distributions in energy and provide fitting formulae for the 
kernels of Faraday rotation and conversion. The coefficients for arbitrary isotropic distributions may be obtained by a further 
integration over energy. Examples are provided for power law distributions covering a finite energy interval and for thermal 
distributions. In the next sections we compare our results for a thermal distribution to results obtained by these approaches. 



Polarization transfer coefficients for a thermal distribution in the high-frequency limit 
Following IShcherbakov |20Qg| ) we express the thermal Jiittner distribution as 

= A e ^fr7j^2 ^ » normalized so that / 47r 7% / 7 2 - 1 m 3 c 3 F (j)dj = 1 . (42) 



where T is the ratio of the thermal energy to the rest mass energy of the particles. By Trubnikov's method, IShcherbakov 
( 20081 ) finds the following polarization transfer coefficients in the high-frequency limit 

h _ Ocosac^jfotr- 1 ) _ f2 2 sin 2 a^ r (K, (T' 1 ) + 6TK? (T' 1 )) 

/hp ~ cw*K,(T-i) ' and hHF_ 2c^ 2 (T-i) ' (43) 

where Ki and K2 are modified Bessel functions of the second kind. In this high-frequency limit his results are identical to 
ours. To compare them, we evaluate / and h from equations (|38|l - (|39[) , also valid in this limit, considering the case of an 
electronic population, so that s q — —1. The / coefficient, derived from equations (|33[) and (55) , becomes in this case 
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log(-h) 




Figure 2. Left: The evolution with temperature T of the ratio of h th to its high-frequency approximation h^ F . Right: The evolution 
of log \h\ versus logT. The T~ 5 / 3 asymptote at large T is a direct consequence of the 7 4 / 3 asymptotic behaviour of the isotropic kernel 
of the Faraday coefficient h at large 7. 



/hf = — T 2 K 2 (T- 1 ) / d7 6 ( 7 £(7) - ^ / ^ T ) ' (44) 



Keeping the integrand part in \Ay 2 — 1 untouched, the term 7 £(7) exp(— 7/T) is integrated by parts, giving successively 

roc roo POO roo 

/ dj 7 £(7) e" 7/T = T 2 / rf 7 - g +T / rf 7 7 p~ 7/T = T 2 ^..(T- 1 ) + / d 7 Tt^T e" 7/T , (45) 

ii Ji V7 2 -l A V7 2 -1 -/i 

where the last equality has been obtained again by integrating by parts the second integral of the middle part. Substituting 
this into equation (|44p , we recover f th in equation (|43[) . Starting from equations (|33l) and (I39|) for /igp we similarly recover 
/ihf in equation (O by noting that / £(7) exp(-7/T) = TAofT -1 ) and that 

/CO 
d 7 7 (2 7 2 - 3)v/7 2 - 1 e _7/T = T {~K 2 {T- 1 ) + & T K3.iT- 1 )) = -TA'o(T _1 ) + 4T 2 (Ki(T~ x ) + 6T ^(T" 1 )) . (46) 



Polarization transfer coefficients for a thermal distribution in the medium-frequency regime 

Stepping out of the domain of validity of the high-frequency approximation, Figure [T] comparatively illustrates the behaviour 
of different approximations to the polarization coefficients by representing their high frequency approximations ()43f) and 
co ntrasting them to th e numerical estimation of their exact counterparts (equations (|25[) and ((26])) and to the fits provided 
by Shcherbakov ( 20081 ) to describe this regime. The difference between the approximations to f th being small, the left panel 



of Figure [T] represents their ratio to their HF approximation. These ratios are displayed as a function of the parameter 

X — T (lO 3 V2 sin a/it) 1/2 , (47) 

defined in Shcherbakov ( 20081 ) . which is proportional to the square root of the left hand side of the inequality in equation 



(|27[) . The right panel of Figure \T\ shows the corresponding approximations to h as a function of temperature, the quantity 
represented in this case being the coefficient itself. This figure uses the following fiducial parameters m = c = 1, a = 7r/4, 
u = 15. It is found that the region-cut approximation to the intermediate regime discussed in Section \3. II and represented by 
the expressions in equations (|31f) and (|32[) provides a good extension of these coefficients to higher temperatures where the 
high-frequency limit becomes insufficient. 



Polarization transfer coefficients for a thermal distribution in the low-frequency limit 

As an illustration of equation (|41[) let us consider again the Isothermal distribution, equation (|42[) . It follows that when most 
particles interact with the radiation in the low-frequency regime we expect 

which corresponds indeed to the asymptotic behaviour that can be computed numerically (using the Olver uniform expansion, 
see section [631 and Appendix [Gj , as shown on figure [2] 
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n= 6- h( 7min ) n=6. 




Figure 3. Polarization transfer coefficients (ieft; fP and right: h? ) for an isotropic power law distribution function as a function of the 
low energy cutoff 7 m i n for different values of the power law index n, as labeled. The coefficient / is normalized to —ns q cosa(ojp r |£2| /cui 2 ) 
and h to cjp r n 2 /caj 3 . This figure has been produced by numerically integrating equations (125l t — (1261) over vj and a for the distribution 
{49}. 



Polarization transfer coefficients for an isotropic power law distribution function 

Let us also consider a power law distribution scaling like 7™ with a low-energy cutoff 7 m , such as: 

Fa (~) = QH(7-7m) (49) 

m7) jV( 7m ,n) 7 "' 1 ' 

For an isotropic distribution, assuming n > 3 for convergence, the normalization factor iV(7 m ,ra) results from equation (|42|l 

V^l7d7 =1 _ (5Q) 



^(7m,n)7 7m 7™ 

The Faraday coefficients can be expressed for such distribution functions in terms of a 7-dependent kernel as in equations 
(|33|l . The low-energy cutoff introduces a Dirac-type singularity in the derivative dFo/d'y at j m , that should be taken care of 
when equations (|33l) are used. Alternatively, these relations may be integrated by parts, the integrated term vanishing because 
Fo vanishes at infinity and the kernels do at 7 = 1. 

The figure [3] shows the polarization transfer coefficients for power-law distributions of various exponents as a function of 
the low -energy cutoff "7111111 ■ These coefficients have been numerically calculated from equations 11251) (126|) by integrating over 
tjj and a for the distribution (|49|l . using equation (|A1|) to express 7 in terms of and w. 

Approximations to the one- variable kernels could have been used instead. We derive in section 16.51 approximations to 
i ?lso (7) and i/ lso (7) obtained by interpolating between their high-frequency expressions and their asymptotic low- frequency 
limits. They provide results of reasonable, though limited, accuracy, but conveniently reduce the calculation to a simple 
one-variable quadrature. The figure [LT1 in appendix [L] compares the results in figure[3]with those obtained from these approx- 
imations. Equations (|33[) . (|50[) and (|40[) indicate that / pl scales for large 7 m as ln(7 m )/7j^. 

In a distribution that is extended in energy, particles of Lorentz factor 7 less than the upper limit y v defined by equation 
(|27l) contribute in the high-frequency regime, as in equations (|38[) and (|39[) , while particles of a larger Lorentz factor bring 
a quasi- resonant contribution that should be added to their non-resonant one. For 7 2> 7 qr , these two contributions add to 
form the asymptotic low- frequency contribution to the coefficients in equations (|40[) and (|41[1 . The QR contribution to the 
coefficient / is small but does not vanish, being the integral of an almost odd function over an interval almost symmetrical 
with respect to zero (appendix [J]). The isotropic low- frequency kernel that results for / has been derived in section [6] and is 
shown in equation (|130|) . 



3.3 Polarization transfer coefficients for anisotropic distribution functions 



The high-frequency results in eq uations (128|) and (J29J) make no assumption about the isotropy, or otherwise, of the distribution 
function Fo- But for the work of Melrosen 1997bh on anisotropic ther mal plasmas, there does not se em to be any other result 
published so far in the literature on general anisotropic distributions. Huang fc Shcherbakov 1 201ll ) give results and fits for a 
number of isotropic distribution functions, including monoenergetic ones. The form in which our results are obtained allows 
to calculate transfer coefficients in the high-frequency approximation not only for isotropic distribution functions but also for 
a large class of anisotropic ones, providing the first terms of an expansion of the anisotropy in multipoles. Strong anisotropies 
still demand a two-variable quadrature, as exemplified by the beam model below. Stepping out of the high-frequency limit in 
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this case would provide one-variable quadratures that are no simpler than the two-variable ones from which they originate. 
We therefore restrict in this section to high-frequency results for different anisotropic distribution functions. 

3.3.1 Quadrupolar and higher multipolar anisotropics at high frequency 

Let us parametrize quadrupolar-type anisotropic distribution functions as the product of a function of the Lorentz factor 7 of 
the particles by a second order Legendre function of the cosine of the particle's pitch angle #, defined such that pu — pcos$: 

F(7,0) = F 2 ( 7 )P 2 (costf). (51) 

For such distributions it is best to switch to the spherical coordinates 7 and The a- and tu-derivatives of the distribution 
function in equations (|25[) - (|26[) or in equations (|28[) - (|29|1 always separate into the sum of two terms, one of which is proportional 
to di ? 2/d7 and the other is proportional to F2. Generically, owing to the factorization property of the distribution function 
as in equation (|51[) . the transfer coefficient could be written, for example in the case of the Faraday rotation coefficient /, as 

/5 n F is ° = Ajf^dn (j± D' f ( % „) + F 2 ( 1 )D f (f, A.)) , (52) 

where Af is a constant factor, still dependent on lj and a, and the functions respectively multiplying ^2(7) an d ^2(7), P/(7 ; /i) 
and Df(-y, fj,), depend on 7 and on the parameter fi — cos a, a being the propagation angle of the radiation (figure[S|. Similar 
quantities refering to the coefficient h can be defined. From equations (|28ll . (|29l) . using equations ()36[) and I) Alfl . it follows 
after integration over $ that, for the quadrupolar model in equation (|51[) 

*ani S o = ns "pr M m 3 c 3 L^^l ( ^(7)(3(2^+ 3 ) P 3 (rt- (-r 2 -1» { 4((^-l)„-(ll^+4)p3 W ) \ 

9 clo 2 J dry \ s( 7 2 -i) 5", v / t 2 -i / 
+ nSq <4lB. mVjyd^i-f) ^(A^ 4 -i>-(*^^ , (53) 

and 

,aniso ^pr^ 2 3 3 f , „, -. / C(-y) ( (24 7 2 +5) P 2 (n) - 12 (2-, 2 + 1) P 4 (nH?) 7< -6 7 4 + n-, 2 +4) +5 (187 4 + 65 7 2 +4 ) (fl) _ 1 2 ( 4t 4 + 37 t 2 + 4 ) p (fl) \ 

^ a>p r fl 2 3^1", ' £ dF% / C( 7 )(l4 7 2 -1087 2 P4(M) + (94-y 2 -7)P 2 ( > x) + 7) 7(4-y 4 -8 7 2 + l<>) + 5( -20 7 4 + 76^ 2 +3l) P 2 (m)+36(2 7 4 -9 7 2 -8) P 4 Q*j N /g^N 
CW 3 / d7 \ 42t( 7 2_!) 2101 / 7 2_1 / ' 

Equations (|53l) and (|54|l are analytical expressions, valid in the high-frequency limit, for the Faraday coefficients / and h for 
an anisotropic distribution function of the form (|51[) . Appendix U gives the corresponding coefficients D'f and Df for higher 
multipoles of the form F n (7)P n (^t) for n = 0, 1 • • • 6. Any axisymmetric distribution F(7,$) can be represented as a linear 
combination of such functions, as F(-y,i}) = F n ('y)P n (cos'&). Since the conductivity depends linearly on the unperturbed 
distribution function, the conductivity of a sum of functions is the sum of the corresponding conductivities. In each individual 
term, neither the multipole "distribution function" Fn,(7) nor the angular factor P„(cos$) need to be positive. The only 
constraint is that the sum over n be positive for any 7 and i?. 

As expected, the coefficients / and h depend on the direction of propagation a in a different manner than those associated 
to an isotropic distribution function given in equations (|38[) -(|39 [) . For a quadrupolar distribution, / varies with a as a 
combination of Pi (cos a) and P3(cosa) while h varies as a combination of P^cosa) and Pt(cos«). 

3.3.2 Beam model at high frequency 

Let us consider now an anisotropic distribution representing a beam. Since theory is constructed for particles, the unperturbed 
motion of which is ruled only by the static magnetic field in the chosen frame of reference, the beam must be assumed to 
propagate parallel to the magnetic field. Otherwise, for example for winds propagating at an angle to the magnetic field, a 
static convection electric field would also be present in the observer's frame and the calculation in section [2~2l would then have 
to be revisited or the Stokes parameters transformed from the beam proper frame to the observer's frame. We form the beam 
distribution F* by boosting an isothermal distribution, so that 

P*(7, 6) oc exp (— 7*(p, 0)/T) , where 7*(p,0) obeys c 2 m 2 ^/l(p, •&) = c 2 m 2 + (p* + p cos(i?)) 2 + p 2 sin 2 (i?) . (55) 
For such a distribution, the derivative of F* w.r.t. a entering equations (|28ll - (l29l ) is given by 

[sj d 2 m? +p 2 + cos(a)p*j . (56) 



dF* = _ F in 



da T uj 7* sin a 

For the purpose of numerical integration, equations (1281) - (1291) and (I56|) are expressed as a function of p± and pu. The result 
of this integration, valid in the high-frequency limit, is shown on Figure [4] for our beam model and for different values of p* 
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Figure 4. Polarization transfer coefficients (/, h) versus the angle of propagation a, represented here on the interval [0, 2ir] for the 
"beam" model of Section 13.3,21 as a function of the beam momentum parameter p* (left panel) and the temperature T (right panel), 
as labeled. The vertical arrows indicate in which sense the parameter varies at ct = 7r. The Faraday response becomes more collimatcd 
along the magnetic field direction as the distribution becomes more anisotropic or colder. 



and T as labeled, with the same fiducial parameters as above. As expected, for p* = we recover the thermal solution of 
Section 13.2.41 but as the beam becomes more anisotropic (i.e. p* increases at fixed T, or T decrases at fixed p* ^ 0), the 
Faraday coefficients are more focussed near the axis of symmetry. This may be understood from the fact that the values of 
the kernels depending on vj and a in equations l]28ft (|29[) are the largest where the difference a 2 — x 2 is the smallest, that is 
where the pitch angle i? of the particle is closest to the direction of propagation a of the radiation, causing a larger response 
when the angle a is in the particle's beam. 

The next two sections formally carry out two of the four integrals involved in equation (|24l) in order to write equations ()25|) 
and ()26|) . and investigate the corresponding resonant and non resonant regimes. A detailed analysis of the various frequency 
regimes allows us to write simple one dimensional quadrature. These two sections can be skipped as a first read. 



4 DERIVATION OF THE POLARIZATION TRANSFER COEFFICIENTS 

We now describe our derivation of the polarization transfer coefficients, which amounts to calculating the elements of the 
hermitian and antihermitian parts of the conductivity deduced from equation (|24[) , that is, of the phase integrals in equation 
(1221) . The usual approach to the phase integrals is to integrate over the time delay r by expanding exp(ia;sin0), or any 
similar expression, in discrete Fourier series such as for example exp(ix sin cj>) — ~^2 n J n (x) exp(in<^>) where J n (x) is a Bessel 
function of relative integer order. An integration over times t' earlier than t (that is over positive delay times r = t — t') is 
then performed that yields resonant denominators ui — k\\V\\ — nfi*. The integrand being proportional to exp(+icjr), and it 
being understood that uu really is a complex Laplace variable, the resonant denominators should be regarded as having an 
infinitesimal positive imaginary part, ensuring the convergence of integrals over r at r = +oo. These denominators are then 
to be understood as complex numbers with an infinitesimal positive imaginary part, the real part of which is the Cauchy 



princ ipal value, the imaginary part being — inSo(u — k\\V\\ — nfi*). This gives an exact expression of the conductivity iBekefil 
19661 ) in the form of series, each term of which involves one such resonant denominator. 



However, since any approximation to the sum of these series has, in the present context, to be carried out up to very large 
values of n, of order cj/ | fi* | at least, this representation of the conductivity is not suitable for our purpose. As when deriving 
the synchrotron emission spectrum, it would be preferable to somehow substitute to t he discret e series a representation in 
which the discrete summation is replaced by an integration over a continuous variable. ISazonov (|l969bl ) attempted this by 



changing the discrete sums into integrals, considering the index n of the Bessel functions as being a continuous variable which 
can be identified with the variable a defined in equation (|23|) . This approach is known to be successful when evaluating the 
dissipative part of the conductivity, in which the many resonant Dirac functions form a sum similar to a Riemann one that 
can be approached by the integral of the so-defined interpolating function, F(a) say. But this is an unfortunate approach for 
calculating the non-dissipative part, in which principal values are involved. A continuous approximation to the derivative of 
F(a) would have been preferable in this case, although any a-priori choice of such an interpolation would be tainted with 
arbitrariness. 



4.1 An exact continuous-spectrum type representation of the conductivity 

We improve here over Sazonov's approach by deriving exact expressions for the phase integrals in equation (|22[) as functions 
of the continuous variable a defined in equation (|23l) . This variable will eventually appear as a continuous index of some 



14 J. Heyvaerts, C. Pichon, S. Prunet, J. Thiebaut 



Bessel functions. This will achieve, without any arbitrariness, the desired continuous-spectrum-type representation. iQin et al 



1 20071) have der i ved th is tr ansformation by using the invariance associated to the periodicity of the unperturbed motion 
while iHevvaertg |l970h and iThiebautl (|201oh derived it by exactly turning the familiar discrete series representation into a 
continuous one by means of an integral represe ntation of the pri nciple values involved in the series. To show here how this is 
achieved, we follow a method similar to that of lOin et al] (|2007h . The phase integrals in equation (|22[) are first expressed in 
terms of the variables in equation (|23[) and the delay time r is replaced by a delay angle ip, so that 



exp(— i s q (acj) — x sin </))) 



s q (jy 



dip 

\n7\ 



exp 



(i (ja + s q e) ip — x sin ip\ ) , with 



in, I 



(57) 



To account for the periodicity of the unperturbed motion, the integration range over the angle ip is separated into segments 
of length 27r, and the integration on the n th segment is carried out by changing from ip to w = ip — {s q <p + 2nir). Then, 



P e (<7, £, 



exp(i (etp + s q x sin (p)) 

Tnj 



J dw exp(i^((7 + 8qS)w — X sin(l(J + 8q(f>)^) . 



(58) 



The fact that the periodicity of the motion is fully taken into account in this exact transformation could be considered 
superfluous because, on the one hand, the period is very long in the ultrarelativistic limit and because, on the other hand, 
the integral over ip in equation (|57[l could conceivably be evaluated by the method of stationary phase applied to a unique 
interval of quasi-stationarity. It will however become clear below that the antihermitian part of the conductivity is in fact 
not determined only by the properties of the motion near the quasi-stationary phase, a situation that may be compared to 
that of a classical plasma, in which the cyclotron resonances determine the rates of emission and absorption but where non- 
resonant particles make the most of the contribution to the medium's dispersive properties. From equation (|58|l we calculate 
the combinations of the phase integrals that appear in equation (J24J) , namely 



/ Po 

{P+-P-)/2i 

\ (P+ + P_)/2 



is q x sin 4> 



2^ e 



dw e 



i(<rw — x sin(w-\-s q (/>)) 



/ l 

sin((f> + Sqiv) 

y C0S(<^> + SqW) 



(59) 



The sum over n of the imaginary exponentials is not absolutely convergent. It should be regarded as a distribution rather 
than as a proper function because it enters an integration over the variable a that will have to be performed at a later stage. 
We define three functions Ti, T2, T3 by 



Ti = 1, 



T2 = sin , 



T 3 



cos 



(60) 



The expressions in equation (|59[) are then inserted in equation (|24[) which gives the components of the current. Then gathering 
the terms wich depend on the gyration angle <p an d integrating over it, the components of the conductivity tensor can be 
written as 

2 



/ dp\\ / p±dp± Mij(px,P\\,kx,k\i) , 

W J-oc JO 

where the matrix Mij is, with i a line index and j a column one: 

( v ± (uj—^-+kiiD(f ))Q 33 (a,x) v ± (uj—^- + kuD(f )) Q 32 (a, x) (v ± uj—^ Q 31 (a, x) ~ k ± vj_D(f ) Q 33 (a, x)) \ 
op± op± 0p\\ 

v ± (uj§^ + fc||D(/ )) Q 23 (a, x) wx(w|^- + k\\D(f )) Q 22 (a, a;) (v±oj^- Q 2 i((T,x) - k ± v ± D(f ) Q2 3 (a,x)) 

v± Op± Op\\ 

y v \\^Wz + h D Uo)) Qi 3 (a;x) v\\{u)-^+k\\D(f ))Q 12 (cr,x) (v^u-^ Qu(a,x) - k x V\\D(f ) Q 13 (a,x)) ^ 



(61) 



(62) 



The coefficients Q a b that appear in equation (|62|l . a and b varying between 1 and 3, encompass the results of the integrations 
over the delay angle and over the gyration angle. They depend on the sign s q of the electric charge of the particle species and 
may eventually be written as 



Q ab (a, x) = (~s q ) a+b 27ry^e 2 " r " CT \ Gab(a, x) , with G a 



d(p f 2n dw 



2tv a 2tv 



T a (^)T h {<p-w)e 
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{4> — w) —sin 1 



The parenthesis in the first part of equation (|63|) can be more explicitly written as 



n \ ^ 2iivno 

27r^e 



(63) 



(64) 



The calculation of the elements G a b involves integrations over the angles <j) and w. lQin et al] (|2007l ) perform these integrations 
along the lines described in Appendix [B] The elements G a b may finally be expressed in terms of Bessel functions of the first 
kind, with indices a or —a and argument x, as listed below: 
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Gn(a,x) = e^JaWJ-afa) G 12 (a,x) = - j -j^(J a (x) J-*(x)) , 

Gi 3 (^) = e'"(-^ + f^)J-»W) G 22 (a,x) = ^e^{^ J' a {x) JL a (x) + , 

(65) 

G 23 (a,x) = g e^A (J CT ( x )J_ CT (x)) , G 3 3(a,x) = sinavr (l - ^ J^x) J_ CT (x)) , 

G21 = — G12 , G31 = + G13 , G32 = — G23 • 

The elements Q a b, and then the matrix elements Mij, axe given in terms of these coefficients in equation (|63|) . The variable 
a being positive, it is useful to eliminate the Bessel functions of negative indices for other Bessel functions with positive 
indices. This is possible since the Bessel function of the first kind and negative index J_ CT (x) can be expressed in terms of 
Bessel functions of the first a nd second kind with positive indices, Jo- (a;) and N a (x) by the left identity in equation (|66[) below 



l Abramowitz fc Ste aun 1964). When derivati ves are involved, the second relation in equation (|66[) . that gives the wronskian 



of the functions of the first and second kind ( Abramowitz fc Steeunl 19641 ). may be used to eliminate N' a if needed. We then 



base the transformation to positive a values on the relations: 

2 

J-a(x) — COSCT7T Jrr(x) — Sm CT7T A CT (x) , J a (x) N' a (x) — j' a (x) N a (x) = . (66) 

nx 

Using equations (|66[1 . the elements Q a b in equation (|63[) can all be expressed in terms of J a (x), N a (x) and their derivatives 
w.r.t. x. Some components of these expressions turn out to be regular at integer values of a, while some other keep a singular 
denominator sin<T7r which generates expressions that must, again, be understood in the sense of distributions. The matrix 
elements My in equations H61|) - (|62p can then be explicitly written as 

Mxx = (+i) lim fl-e M ") ^ + k „£>(/„)) n - 2 ( 4(x) - Mx)N„ (x) - — ) 

MxY = (~s q ) lim ( 1 - e** 1 ") ^-r + k»D(f )) ■K-(^^J t r{x)j' a {x)-j' tT {x)N tT {x)- 

JV->oo V / | 12* | Op± X \ Sill (TIT 
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Spy 

Let us now specify the meaning of lim (l — e 2l7riv<T ) , which enters with other factors in integrals over a implied by the 
integration over momenta in equation (|61|l . When this factor multiplies a regular function of a, its limit is unity since, when 
N diverges, exp(2i7rA r cr) oscillates infinitely rapidly leaving in the limit a vanishing integral when multiplied by any regular 
function. The case when a factor (cos a/ sin an) multiplies (1 — exp(2inNa)) deserves closer scrutiny, since both sin<T7r and 
(1 — exp(2i7rAo-)) vanish at integer values. It is shown in Appendix [C] that the limit of their product is the distribution 



N 



Um ( 1 _ e **»A S2^JL s V (S2^1) = P ( c ^l)^ i y SD (a-n) . (69) 
\ / sin an \ sin an / \ \sincr7r/ z — ' / 

where N is the set of positive integers, So a Dirac distribution and the notation V means that the function which follows in 
the parenthesis should be taken as a Cauchy principal value near each of its singularities, which in this case are the integer 
values of a. The effect of the multiple resonances at all integer values of a that are involved in the matrix elements of the 
conductivity are concentrated on those terms in equations (|68|) in which factors (cos an j sin<T7r) subsist, although it could 
have been expected that resonances should be present in all terms since equation (|64|) exhibits such a singular factor. In 
changing according to equation (|66[) for Bessel functions with positive indices only, some of these singular factors have been 
disposed of, having been regularized owing to the favourably phased term proportional to sin an N a (x) in equation (|66p . This 
means that the transformation to Bessel functions with only positive indices has effected the summation over resonances for 
these favourably phased terms. Resonances only remain explicitely present in the unfavourably phased ones, that originate in 
the term cos an J a (x) in equation (|66[) . It is shown in section l4~4l that these residual resonant terms eventually turn out to be 
negligible when a;/|f2|^> 1. 
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Figure 5. The Z axis is along the unperturbed magnetic field. The z' axis is along the wavevector k = fccoscte^ + k sin a ex, where 
< a < 7r. The axes Y and x' arc identical 



4.2 Polarization transfer coefficients from the conductivity 



The elements of the transfer matrix are given in equation 1)121) in terms of the components of the conductivity tensor in the 
plane perpendicular to the direction of propagation of the radiation. Since the components of this tensor are known in a frame 
in which the static magnetic field is along the Z axis, the components in equations ()68[) must be transformed to the new 
reference frame x' , y' , z' represented in Fig. [5] in which the wave vector k is along the z axi s. For our purpose, it suffices to 
calculate the transverse components x'x' , x'y' , y'x' , y'y' which are given bv lThiebautJ l|2010h 



M,, T / = My 



M x 'y' = sm a Myz — cos a Myx 



M y i x i — sin a Mzy — cos a Mxy 



My>y> 



cos a Mxx + sin a , Mzz — sin a cos a (Mzx + Mx2 



(70) 
(71) 



These calculations are straightforward, but it nevertheless takes some algebra to reduce them to the simple form in equations 
(|72[) (|74[) . Some details are given in Appendix [D] The transverse components of the M matrix can eventually be written as: 
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(72) 
(73) 
(74) 



From equations (|72 )l -(|74 p , we may calculate the corresponding components of the conductivity. The density n of the particles 
is factored out of the distribution function by writing it as /o = nfb, defining a reduced distribution function normalized to 
unity and a plasma frequency uj pr for the relativistic species of particles considered: 

„2 
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2np± dp_ 
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nq 



(75) 



The components of the antihermitian and hermitian parts of the conductivity can then be written as 
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a$ y , 2 a&fi 2 ^ ffm 3 c 3 2 8F - 

= -27T p ,, > // — 5— dwda zu — — JJx)So[a -n) . (81) 
£0 w 3 sm'a 9cr 

The elements of the transfer matrix in equation (JTJ can be found from these results by using equation (|12p . Our choice of 
reference axes in the plane perpendicular to k results in the vanishing of the coefficients Kqv and K/y, since a^, y , is real 
and cr^/yi is imaginary. In equation (|78[) the integration over zu and <r of the terms cos a(crd TJ7 i ? o + zudaFo) can be reduced to 
an integral on the boundary of the physical domain which vanishes. 



4.3 Formal expression of the polarization transfer coefficients 



The components of the hermitian part of the conductivity in equations (|79|1 - (|81[) are associated with dissipative radiative 
effects. Equations (|12p relate them to the absorption coefficients per unit time which appear in equation |l}. The transfer 
coefficients per unit length in a medium of negligible dispersion are obtained from them by dividing by the velocity of light. 
The discrete sum over the large positive integers n in equations (|79|l (|81[) can be approximated by an integral, so that 
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We return to the dissipative coefficients in section 15721 There are two non-dissipative coefficients, which in a stationary medium 
are usually defined per unit length. Assuming again no dispersion, they are / = Kqu/c and h = Kuv/c. When acting alone 
on radiation propagating in the direction of the unit vecor n, they cause the Stokes parameters Q, U, V to vary according to 
the equation: 
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(85) 



This simplified transfer equation leaves the intensity I invariant as well as the degree of polarisation \/ Q 2 + U 2 + V 2 /I. The 
argument of all Bessel functions implicitely being x, the coefficients / and h can be written as 
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(86) 
(87) 



4.4 The residual contribution of resonances to non-dissipative coefficients is negligible 



The principal value terms in equations (|86|l (|87|) are negligible in the limit uj 3>| S7 ] in which a and x are large. To prove this, 
we use the fact that the interval between successive zeroes of sincr7r is unity, whereas the Bessel functions vary on a much 
longer scale because both their index and argument are large, the argument remaining smaller than their index, though. This 
can be seen from the definition of a and x in equation (|23l) . The contribution to the integral over a of each unit interval 
[n — ^,n+ |] can then be calculated at any accuracy by Taylor expanding about n the function which factors cot<T7r in 
equations (|86[) and (|87[) . This provides the result of the integration on [n — \, n + J] in the form of a series. The summation 
of these functions of n over all unit intervals centered on integer values can then be approximately replaced, provided a is 
large, by an integral over a since they slowly vary with n. This integral turns out to be extremely small, owing to particular 
properties of the functions J a (x) and J'&(x) for large a and small x. Details are to be found in Appendix[E] With this further 
simplification, the polarization transfer coefficients in equations (I86|l and (|87[) lose their principal value terms. Terms in which 
no Bessel functions are involved can be reduced to a line integral on the boundary B of the physical domain in the zu-a plane 
(Appendix [A} , which is particularly useful to transform equation (|87p in which the term devoid of Bessel functions is the 
divergence of the vector with components = zuFo and V a = —aFo. All calculations done, we get 
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The notation ai,(zu) = \fa^ + to 2 denotes the value of a at a point of abcissa vj on the boundary B of the physical domain in 
the w-a plane. It will soon be shown that, when expanding the Bessel functions in \Q,\ /uj, the non-Bessel terms in equations 
(|88l) and (J89J) almost exactly cancel the zeroth-order terms of the development. Equations (|88|) and (|89l) are applicable to 
any distribution function. The integrals over a and/or zu extend to the full range of physically relevant values and no other 
approximation has been made than neglecting the series of residual principal value terms, which is well justified. Equations 
(|88|l and (|89|l are then close to being exact. They however feature complicated kernels, which makes it desirable to find simpler 
and suitable approximations to them. 



5 NON-RESONANT AND QUASI-RESONANT PARTS OF THE TRANSFER COEFFICIENTS 
5.1 Quasi resonant contribution to transfer coefficients from Nicholson's approximation 

The phase integrals in equation (122C involve the integration of trigonometric functions which may vary more or less rapidly 
with the delay time r, depending on the values of the particle's parameters and the gyration angle <j). Let us denote by i? the 
pich angle of a particle, that is, the angle of the particle's velocity with the magnetic field. According to equation (|22[1 . the 
characteristic variation frequency of the phase most often is of the order of uj. When however the modulus of the particle's 
velocity is close to the speed of light, the characteristic frequency of phase variations may occasionally be much less, when the 
angle between the particle's velocity and the wave vector becomes small enough. For example, equation (|22[) indicates that 
this frequency would be of order | Q | or less when 

(V \ V V 1 11 

1 M — (1 — cos(i9 — a)) H — sin a sini? (1 — cos <j>) < - — -. (90) 
c/ c c ck 

Since u = ck/\ f2| is very large, this requires that each term on the left of the second inequality in equation (|90|) be less than 
the term on the right, i.e. that 

i«-*i<(m"\ w<(m"°, (w < nay. m 



For quasi-resonance, the third inequality in equation (|91[) requires that uj should be less than 7 2 \Q\ which is of order of the 
characteristic frequency of the synchrotron emission spectrum by a particle of Lorentz factor 7 (equation (|101[0 . The other 
inequalities in equation (|91|) require that the angle between the wave vector and the particle's velocity be, for a frequency at 
the peak of synchrotron emission, less than I/7. These inequalities are only satisfied by a restricted class of particles, which 
we refer to as quasi-resonant particles. For quasi-resonant particles, the characteristic frequency of the variations of the phase 
(the integrand in equation (J22J) i s i during a small fraction of the synchrotron gyration period, much less than uj and, because 

of this slow variation with r, the phase integral is exceptionnally large. 

This induced some authors |Sazonovlll969bl ; lHevvaertslll970l ; rMelroselll997c1 ldh to consider that quasi-resonant particles 



entirely det ermine the no n-dissipative transfer coefficients, just as they determine the dissipative ones and the emission 
coefficients |WestfoIdlll9"59T l. This shows up when phase integrals are evaluated by the method of the (quasi-) stationary phase, 



resulting in the presence in the results of Airy-type functions. The dominance of the contribution of quasi-resonant particles 
to non-dissipative transfer coefficients is however not granted. The functions which determine the dissipative coefficients in 
equations (|79|) (I81|l happen to be vanishingly small for non-resonant particles, but those determining the non-dissipative ones 
have a much wider support and are associated with a large number of non-resonant particles. At this point, we do not know 
whether the quasi-resonant contribution to the non-dissipative coefficients is negliglible compared to the non-resonant one, or 
comparable to it. This will be discussed in section [5741 

The first inequality in equation (|90p . which defines quasi-resonant particles, can be written for (f> = in terms of the 
variables a and x in equation (|23p as (a — x) < 7. The pitch angle # of the particle and the propagation angle a of the 
wave being almost equal at quasi-resonance, a ~ 711 sin 2 a. Considering the general case when a is not very small, this yields 
a ~ 711, a very large value. Imposing the condition (a — x) < 7 to a wave with frequency in the peak of the the synchrotron 
emission, such that uj ~ 7 2 fi|, we get a ~ 7 3 , which then implies that for quasi-resonance 

Thus, for quasi-resonant particles, a and x are very close to each other, x being however slightly smaller, since, from equa- 
tion (|23l) . (a — x) must be positive. The inequality in equation (|92|) places the argument x in the inter mediate r egion where 
th e so-ca l led N icholson's approximation to Bessel functions of large index and argument is appropriate ( Watson ( 19221 ') chap. 



Olverl (|l952r i). The Nicholson's approximation to J a (x) is well known but the corresponding approximation to N a (x) is 



8, 

not. Its derivation is outlined in Appendix [F] The functions J CT , N a and their derivatives may be represented in this range 
by modified Bessel functions of the second kind K v (g), with index v = 1/3 or 2/3, and by a similar combination L v (g) of 
modified Bessel functions of the first kind, I± v {g). K v and L v are defined by 
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Figure 6. The physical domain of the ro-cr plane is above the hyperbolic line £?, of equation <r 2 = to 2 + oq. The smallest possible value 
of a is (Jo = uj sin a/ \ Q \ . The quasi-resonant domain is the shaded area defined by g < 1 (equation d 94 1 ) ) . It is bounded by the dotted 
line Bqr, represented by equation 1971 . The smallest value of a on £3qr, <r qr , is given by equation I I97H for to = 0. The oblique line is 
the locus of a constant value of the Lorentz factor 7. It intersects B at ro = ro_(7) and ro+(7) and £?qr at to = ro gr _(7) and ro 9r +(7). 



2 sin i^7r 



(93) 



The definition of the functions K v is standard (jAbramowitz fc Steeunlll964l ) and applies to integer values of the index in a 
limit sense. Our definition of the functions L v in equation (|93[) does not make sense for integer v. However, we only deal here 
with indices v — 1/3 or 2/3. The argument g on which these functions depend is 



From equation (|92l) . it can be seen that the variable g is 0(1) or smaller wherever the Nicholson's approximation applies. The 
latter yields the following approximate representations: 



Ja{x) 
J'„(x) 
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a — x 



Ki(g), 



(95) 
(96) 



The quasi-resonant contributions to the Faraday coefficients / and h are obtained by substituting equations (|95H (I96[) into 
equations (|88ll H89p . the domain of integration in the uj-o plane then being restricted to the quasi-resonant one. This domain 
is characterized by the inequality in equation (|92[) . or equivalently, since at quasi-resonance x and a are nearly equal, by the 
requirement that the variable g in equation (|94[) be less than or equal to unity. From equations (|94[) and (|A2|) . the condition 
that g = 1 translates into 



2 . 2 Q 2/3 4/3 



(97) 



The variable g can then be less than unity only when <r exceeds a threshold cr qr and the Lorentz factor exceeds a related one, 
7q r , such that the line 7 = 7 qr in figure [6] be tangent to the boundary Bqr. The Lorentz factor 7 qr is of order of the 7 variable 
associated with cr qr and zu = 0. To sum up, 



3/2 
£o_ 

y/3 



7qr 



3 sin a 



(98) 



Thus, a particle can only quasi-resonantly interact with the radiation when <r exceeds <r qr . Since oq defined in equation 
is usually large, <r qr ^> <to. The characteristic frequency w c (7, a) of synchrotron emission by particles of Lorentz factor 7 qr , 
defined in equation (|101[) . is almost w. The figure [5] represents the quasi-resonant domain in the w—a plane. When u — uj/ \Q,\ 
becomes very large, <r qr diverges as it 3 / 2 while the <r value corresponding to a given Lorentz factor 7 diverges as u. This means 
that at high enough frequency there will be a negligible number of particles in quasi-resonance and the transfer coefficients 
will then be essentially given by the non-resonant contribution. When on the contrary a > cr qr for relevant values of 7, the 
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quasi-resonant region partly contributes to / and h. These contributions are: 

u 2 Q. 2 ff m 3 c 3 , , f 2V2 f(T-X\i 7T A OFq 

te = ^±rJU*K da '* rio *{ife (— ) (99) 

h QR =^[[ Hffdu,da( 4 4 (^) 2 Kl (g)L,( g )+ 2 -^ (^)K l{9 )L l{9 ))^. (100) 

cw 3 JJ QR sm 2 a V\/3 V x ) 3 3 ^ ^3 V x ) da 

where the subscript qr indicates that the integration should be carried over the quasi-resonant domain only. 



5.2 Polarization-dependent absorption coefficients 

The dissipative coefficients per unit length are given by equations (1821) (I84[) . Since the product of Bessel functions declines 
exponentially out of the quasi-resonant domain, it is appropriate to use the Nicholson's approximation to represent them 
(equations (|94[) . (|95[) . (|96[) ). Equations 1)821 1 1 )84)1 are then expressed in terms of the Lorentz factor 7 and the pitch angle 1? of 
the particles, or the difference ip of the latter to the propagation angle a, which remains small in the quasi-resonant domain. 
This results in d a Fo « (d-jFo — ip d#Fo/y)/u sin 2 a. Owing to the smallness of ip, the angular derivative term may sometimes 
be neglected. Similarly, the argument g of the A Bessel functions in equations ()95)) and ()96|) may sometimes be approximated 
by So ~ (o>/37 2 fi sin a) (1 + j 2 ip 2 ) 3 ^ 2 . The c ritical frequency to c of synchrotron emission by a particle of Lorentz factor 7 in 
the direction a is defined by (| Westfol J 19591 ) : 



| 7 2 |fi| sin a. (101) 



The in tegration over the particles directions of the right hand sides of equations ()82j) 1)84) 1 can be performed as in IWestfoldl 
( 19591 ). which gives for the total and linear polarization absorption coefficients: 



Ku W 2 r |n| y/i-K . f 3 3 2 , dF LU 



sm a m c 



-rd-y-fT- — K 5/3 (u)du, (102) 



c ecu 2 2 J dj lj, 

KlQ Wpr 1^1 ^3?r . f 3 3 2 , Wo w / OJ , 

' mc -y d-y — — — A 2/3 — . (103) 

O"/ U)c I Uc ' 



Here, d ~,Fp is meant to be taken at ■& = a, the propagation angle. The results ()102|) and ()103|) coincide with those of ISazonov 



(1969a), considering his use of the CGS system of units, the definition of his distribution function, and the presence of an 
unfortunate typo in the first of his equations (2.2), where f^J v A s / 3 (u) du should be replaced by jf- / , A 5 / 3 (m) du. 

The calculation of the absorption coefficient for circular polarization in equation (|84)l is less straightforward because its 
dominant order contribution involves the integral of an odd function of zu, or ip, which vanishes. Thus Kiv/c generically is 
much smaller than the other two absorption coefficients. Using the Nicholson's approximation for J a (x) and J' a (x), equation 
(|84[l becomes: 

K IV 4V2 u* r n 2 ff a 33. — dF (a - x) 3 ' 2 

= 5- s q 3 u m c smi9V7 2 -ldydip — — vj ± K 1/3 (g)K 2/3 (g) . (104) 

c 3 cu A J J da x 1 -' 2 

The integrand on the right hand side of equation (|104)l should be expressed in terms of 7 and $ from equations ()A1[I and 
expanded to the first non- vanishing even order in ip. This implies that the angle derivative term in d^Fo (equation (IA6)) ). 
which is of order ip, should be accounted for, that the derivative O^Fq be Taylor expanded about 1? = a and that all other 
factors involved in equation ()104)) be similarly expanded. This applies in particular to the product of Bessel functions, owing to 
the fact that the actual value of their argument g, given by equation (|94)) , slightly differs from its lowest order approximation 
in ip. With two terms, this argument g may actually be expanded as 

9 = 9 o + * = ^ ( ^ + 4) 3/2 - * f^ 2 + i) ^ • ( 105 ) 

3sma \ 7 i y 6 sin a \ TV 

Integrating the correction to the product of Bessel functions over ip by parts, these expansions lead to 

Kiv 2 o, 2 r ffm 3 c 3 3 <9A cosa / 2 2 W 2 l\ 3/ V /■ Mr f 



3 "c 



where, again, all deri vatives of Ap are to be taken at $ = a. The integration over ip is carried out following a procedure similar 
to that described by Westfold ( 1959T ). which gives 



73 
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*q- 



K 1V 2tt u> 2 w [Oj /■ S3 6>A / I"" , . w „ (u 
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sina|m 3 C 3 7 d 7 (J^ - ) 



(107) 



The result in equation (|107|) coincides with that given bv lSazonovl (|1969al ) in his equation (2.2) 



5.3 Non-resonant contribution to transfer coefficients from Debye's expansion 



In the non-resonant regime, the reverse inequality to equation (|97[) applies, though x and a both remain large but need not be 
almost equal. These conditions a re suitable for using the Debye expansion of Bessel functions of large indices and argument 
( Watsorj|l922l ; lMatvivenkolll992l ). subject to the condition that (a — x) ^> <t 1//3 and, in our case, that x < cr. The Debye 
expansion represents the Bessel functions of large index a for a given value of the argument-to-index ratio x/a, represented 
by a parameter £ such that 

- = -^. (108) 
a cosh £ 

The Debye expansions of J a (x) and N <T (x) can be written as |Watsonlll92i ): 

J a (x) - J (y ^ cqsH J - + V 27r(jtanhc r(1/2) ((7t anhC)™' ( ' 

N(X ) = N ( ° \ V2e^-^ ^ r(m + l/2) (-ir2- A m (Q (uo) 
ctW " HcoshJ V27ratanh^ ^ T(l/2) (<rtanh£)™ ' liiUJ 

where the coefficients A m depend on f, but for Ao that equals unity. T(y) denotes the gamma function of argument y. The 
n-th term in the sum is of order a~ n , that is of order (|fi| /ui) n . Here we only need to proceed to second order. It can be 
checked that at this order the Debye approximation continuously merges into the Nicholson's approximation at their common 
limit of validity. While a is regarded as a large parameter in the expansions in equations (|109|l and (|110|1 . the ratio x/a should 
be considered of order unity. Approximations to J' a (x) and N^(x) are obtained by derivating equations (|109|1 and (| 1 10[l with 
respect to x at fixed a, taking care of the fact that the auxiliary variable £ depends on x, and thus also the coefficients Ai 
and Ai- Some terms resulting from this derivation contribute at this order. All calculations done, we obtain 

Mx)N a (x) = . ( 1 + , (HI) 



7TVO 

f( v Va 2 -a- 2 / 6^2+2^ 2Ai x 2 A\ 1 x 4 \ 

j CTW iv aW - ^ ^i-r ff2 _ a;a + (ff2 _ 3,2)2 4(a 2 - a: 2 )V- 1 J 
where 

- i _ A ^ 4 JL _ c^ 2 385 a 4 , _ __5_ aV 

1 ~~ 8 ~ 24 a 2 - x 2 ' 2 ~ 128 ~ 576 cr 2 - x 2 + 3456 (a 2 - x 2 ) 2 ' 1 ~ 12 {a 2 - x 2 ) 2 ' { ' 

The first terms in the parentheses of equations (|111[) - (|113|1 . that are equal to unity, remain in the unmagnetized limit, — ^ 0. 
The non-Bessel terms in equations H861) - (|87|) are com P ara ble to them. The non-resonant contributions to equations 1)88[) and 
(18911 are obtained by integrating these expansions of the Bessel functions over the non-resonant domain, denoted by the suffix 
nr: 

, 2tt% o; 2 r fi 2 ff m 3 c :i w dF ( 1 x 2 6A 2 + xA[ - A\ 3 A x x 2 \ 

/nr - — —JJ m — ^ drod<T - -fa [2 (a 2 -x 2 r/ 2 + a^x~ 2 + 2 (a 2 - x 2 ) 2 ) ' (U5) 



n 2 u 2 D. 2 ff m 3 c 3 , , dF 2Va 2 - x 2 A 6^2 — A? + zAi Aiz 2 1 

/ . 9 drorfa — ( 1 H h 

/nr 



coj 3 JJnr sin 2 a 9cr 7r \ cr 2 — a; 2 (a 2 — x 2 ) 2 8 (cr 2 — a; 2 ) 3 

3„3 ap. „2 / cyl. /|2\ „,,2 n 2 r+»™3J 2 . 2 



^'j w dFo a 2 ( , 6A2 - Af\ ^ 2 r » 2 f + °°m 3 c 3 2 ro 2 + q 2 

3—// ^-d-wda — 1 + — — + — i^— / - T ^-dro— F (tu,tT 6 (CT)). 116) 

JJnr sm a 00 -k^ct 2 - x 2 \ o 2 - x 2 J cu 6 J _ ao sm 2 a ^/vj 2 + a% 

It is remarkable that the terms independent of magnetization have canceled out of equation (|115l) . They also cancel out of 
equation (|116[) as we now show. Gathering and arranging them, these terms can be written as 



(0) 7wg r 2 ff m 3 c 3 9F 2B7 a + gg , nu 2 pl tf f+°°m 3 c 3 j 2m 2 + a 2 

h NR= 3—// ^-^-dvoda — + — / -^-j-eto — === F (n7,cr b (w)) . (117) 

"^ 3 JJnr sm 2 a da ^ 2 + a J cu; s 7_ oc sin 2 a ^/ro 2 + cr 2 



The double integral term may be integrated explicitly over a at given vo since only 9 CT fb depends on this variable. The 
integration is on the values of a between its value o-&(ro) on the boundary B of the physical domain and its value aQn(vo) — 
(uj 2 + <To) 3//4 /v / 3 on the boundary Bqb, between the non-resonant and quasi-resonant domains. The contribution from the 
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lower boundary, at(zo), cancels the term that was already in the form of a single integral over zu. The contribution from the 
upper boundary at aq^izu') remains and could as well be considered to be part of the quasi-resonant contribution since it can 
be written as an integral over the quasi-resonant domain. Naming this contribution /ibqr: 

nuj 2 pl .n 2 [+°° m 3 c 3 2zu 2 +a 2 jrw^O 8 ff m 3 c 3 2^ 2 + a 2 9F 

/ibqr = — / -^-^-dzu Fo{w,a bqr {w)) = — // dzuda —=== — . (118) 

cw d J_ oa sma ^/zu 2 + a 2 cw 6 JJ Qn sma y/zu 2 + a$ oa 

We refer to the remainder of the non-resonant contribution to h as the reduced non-resonant contribution /inr: 

- cupfl 2 r r m 3 c 3 f n f 6A 2 -Aj + xA' 1 A x x 2 1 x 4 \ 2 ( &A 2 - A 2 \\dF . 

kNK ~ 77 sin^ dwda K \ (* 2 ~ * 2 ) 1/2 + (* 2 ~ * 2 ) 3/2 ~ 8 (<r» - ) ~ a ° [ {a 2 _ x 2 )3/2 ))-g^- ( 119 ) 



The Sqr contribution to h in equation (|118[) is then associated to the quasi-resonant contribution /iqr in equation (| lOOf) to 
form a reduced quasi-resonant contribution /iqr: 

/iqr = ^dzuda (^(^) 2 K 2 m , {g) +*£(fLZ*) Kl (g)L ,(,)-„ ^ + ^ ) f° . (120) 

The combinations of the functions Ai, A 2 and xA\ that appear in equations (|115[) and (|119[) may be expressed in terms of a 
and x by using equation (|114[) . When the integration over the NR domain can be extended to the full domain, this results in 
equations 



5.4 Non-resonant versus quasi-resonant contribution to the Faraday coefficients 



The quasi-resonant contributions to the Faraday coefficients are given by equations (|99[) and (|120[) and the non-resonant 
ones by equations (|115[) and (|119|) . We now compare them, making a simple Ansatz concerning the distribution function, 
namely that Fo is a linear function of zu, so that d a Fo = Foo-(cr) + zuFo aTa (a). The Fq& term produces a nil contribution to 
/ if, as assumed here, it depends on a alone. The zjFo am term produces a vanishing contribution to h since its kernel is an 
even functions of zu. This Ansatz is sufficient to give a hint on the dependence of the transfer coefficients on the distribution 
function. The contributions /kr and /qr to /, or /inr and /iqr to h, may then be written in a form involving a kernel 
depending on a, such that 

fn = 2s q ^-f ^-daK^\a)F ^(a), = f ^ da K^(a)F 0a (cr) . (121) 

cu) A j R sin a cu A J R sin a 

where R stands either for the QR or NR domain and (<r) and K ^ (a) are corresponding one- variable kernels, integrated 
over zu at given a. For Faraday rotation: 

- fj™ \ , = jj* % *?* (k V3 ( 9 )L 1/3 ( 9) - . (122) 

The quasi-resonant kernel for /, on the right of equation (|122|l . can be expressed in terms of the variable g from equations 
(|H4[l and then be integrated over the quasi-resonant domain g m {o) < g < 1. The lower bound (/ m , defined in equation (|H4[I . 
approaches zero when a diverges. The integral on the right of equation (|122l) converges as g m approaches zero, which it does 
when a grows much larger than a ql . The asymptotic non-resonant part of the kernel of / can be calculated as shown in 
appendix [H] The results, valid for a cr v , are: 

*£2W « ln (f) . K qI(°) « °- 6 ^ 2 ■ (!23) 

The non-resonant contribution to / slightly dominates when a 3> o" qr - These results partly depend on our assumption that 
the distribution function depends linearly on zu. Otherwise, /nr would also have a term proportional to Fo a that is absent 
from equation (|121l) by imparity. We return to this in section ROl 

The zu- integrated kernel for h being defined as in equation (|121|) . its non-resonant part results from equation (|119l) and 
is calculated in Appendix [H] The quasi-resonant part /iqr is given by equation (|120|l and its calculation is also outlined in 
this Appendix. The kernels -Knr(ct) and Kq^((j) turn out to be approximately given, for a 3> f qr by 

vWt ^~ 77 ( a \ 4/3 isW, v ^ . t 4/3 

^NrW W -g VgJ ' ^Q R (c r )~+2 Cr ■ ( 124 ) 

A glance at equation (| 124[) shows that both contributions to the kernel vary when a 3> o- qr as cr 4//3 , the quasi-resonant 
one being numerically larger. 

The scaling in a 4 ^ 3 anticipated for the h kernel from equations (|124[1 does not c onflict with the fact that th e Faraday 
conversion coefficient h decreases to zero at high temperature for thermal distributions ( Huang fc Shcherbakovll2011 ). We have 



confirmed this decline of h th with temperature by numerically integrating equation (|89[) for a thermal distribution function 
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and checked its compatibility wit a kernel in 7 4,/3 . We return to this is section [6721 The Faraday conversion coefficient h is 
largely affected by the quasi-resonant contribution to it, which is of an opposite sign to the non-resonant one. Equations (|31[l 
and (|32[) propose a simple, though crude, way to account for the increasing importance of the quasi-resonant contribution 
which should be good enough wherever the balance passes from non-resonant to quasi-resonant domination. For a given uj, the 
QR domain becomes narrower in angle as 7 increases and the QR kernel has an effectively compact support in this domain. 
Therefore our calculation in equation (|124|) of the quasi-resonant contribution to h should be close to being exact at large 7's, 
since dFo/da certainly is almost constant with zu over the narrow QR domain, as was assumed in this section. 

The condition for the quasi-resonant domain not to contribute to the transfer coefficients is that the distribution function 
be negligible in it. From equation (|92|) . this happens when for most particles a — \Zu 2 — <Tq > cr 1 ^ 3 , or equivalently when the 
variable g in equation (|94p is larger than unity. Expanding the expression of the latter condition in 00 /& and noting that 
a = 7itsin 2 a at zu = 0, this inequality reduces again to the condition that equation (|27p be satisfied. Up to a somewhat 
arbitrary factor, the paramet er on the left of equ ation (1271) is similar to the s q uare of the regime-change parameters X and 
^qXa respectively defined by IShcherbakov |2008l ) and iHuang fc Shcherbakovl (j201lh . When the inequality in equation (|27|l 
is not satisfied, the coefficients / and h differ from their high-frequency approximations in equations (|145|) and (|146|) below, 
where the integration is meant to extend over the full zu-a domain. The intrusion of quasi -resonant contributi ons is at the 
origin of the change in the trend of the variations of / and h with temperature described by IShcherbakov for thermal 

distributions. This author attributes this change to a failure of the high-frequency expansion, which is correct in the sense 
that the quasi-resonant contribution cannot be expressed in the form of a series expansion, no partial sum of the Debye series 
being able to represent the Bessel function in the domain of validity of the Nicholson's approximation. 



6 LOW-FREQUENCY KERNELS OF FARADAY COEFFICIENTS OF ISOTROPIC DISTRIBUTIONS 
6.1 Remarks on the low- frequency limit of isotropic kernels for the Faraday coefficients 

In the limit of large 7's, the quasi-resonant domain contributes little to the Faraday rotation coefficient /, given by equation 
(|88| ), because it involves the integral of an odd function of the difference ip — # — a on the small QR domain. This remark holds 
also for the integral over this same quasi-resonant domain of the non-resonant contribution. For an isotropic distribution, there 
is however a non- vanishing NR contribution to / that is proportional to d 1 Fo. Such a contribution is absent for distributions, 
as considered in section l5~4l that are linear in zu. 

By contrast, the quasi-resonant domain, if substantially populated with particles, largely contributes the Faraday con- 
version coefficient h. Particles in the QR domain satisfy the inequality inverse to that in equation (|27[) in a strong sense. We 
refer to the case when this strong inequality is satisfied as the low-frequency limit. Since any isotropic distribution function 
depends on 7 only, the Faraday coefficients for such distributions can be reduced to a single quadrature over the Lorentz 
factor involving a 7-dependent kernel as defined in equations (|33[1 . The exact expression of the isotropic kernel for / results 
from equation (|88|l in which the variables zu-a should be changed to zu-^y, giving on integrating over zu at a given 7 

F i 80 = _2tt% ^ x fj^ x)N ^ x) + ±_\ (125) 

sm'a J m _ V "/ 

The exact isotropic kernel for h is similarly obtained from equation (|8T[) . in which the principal value term should be neglected. 
Using equations (|A1[) - (|A7[) . the isotropic kernel for h can be written as: 

H is °h) = -^r- I + dzu(x 2 j' a (x)N' a (x)-zu 2 J a (x)N a (x)--( 1 usii?a + 2mcosa)\ . (126) 
sm a J w _ \ ty ) 

In equations (|125|) and (|126[) . the bounds zu±('y) are defined by equation (|37[) , The kernel H ls °(^) could have been deduced 
as well from equation (|89[) . reformulating the line-integral term in it as a surface integral. These different expressions are 
equivalent because the non-Bessel integrals in them, though apparently different, are actually equal, given the particular 
values of the bounds. 



6.2 Thermal Faraday conversion coefficient in the low-frequency limit 



To check the behaviour of the Faraday conversion coefficient h in the low-frequency limit, we have numerically integrated 
its quasi exact expression in equation (|126f) with Mathematica for a thermal distribution function. A t a giv en frequency, the 
low- frequency limit would then correspond to the limit of high temperatures. Huang fc Shcherbakov 1 201lh have found that 
in this limit the Faraday conversion coefficient decreases with temperature. We confirm this and find that this decrease scales 
as — J 1-5 / 3 by extending our numerical calculations to high values of the temp eratur e T (se e figure f2}. This has been possible 
by using the uniform expansion for high order Bessel functions derived by Olver ( 19541 ) which allows a fast and reliable 
calculation of high order Bessel functions (see Appendix [G) . Returning the result of the integration over the directions of 
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# 1S0 (r)/(y 4/3 sin(a) 2/3 ) 




Figure 7. The evolution of H lso (7)/(7 4 / 3 sin 2 / 3 «) (computed via equation (1 1 26 p using the Olver expansion) as a function of 7 for 
different values of a. The radiation's parameter u = 15. For all values of a, this ratio asymptotes to a constant for large 7. The dashed 
line corresponds to the asymptotic prediction of equation ^ 1381 - 



the particles before performing that over 7 reveals the high energy properties of the isotropic kernels derived in sections 16,31 
and 16.41 This numerically supports the analytical result, reached in section [6.41 and shown in figure [7] that the kernel of h 
asymptotically grows with the Lorentz factor as 7 4//3 and shows that this increase with energy of the kernel is consistent with 
the fast decline with temperature of the thermal Faraday conversion coefficient. Note finally that, considering equations (|42|l 
and (|33[) , equation ([138} predicts at large 7 the result of the numerical integration and the J 1-5 / 3 asymptote for the thermal 
distribution. 



6.3 Isotropic kernel of the Faraday rotation coefficient in the low-frequency limit 



The integral in equation (|125j) may be separated into non-resonant and quasi-resonant parts as follows. The Faraday rotation 
kernel i ?lso (7) is the sum of the expression in equation (|115[) , integrated over the NR domain and limited to its lowest order 
term in Q/uj, and of the QR contribution from equation (|99[) . The integral over the NR domain in equation (|115[) can be 
extended to the full domain at the expense of adding, if needed, to the contribution from equation (|99p a correction to 
account for the undue integration over the quasi-resonant domain so introduced. After the necessary change of variables, and 
considering the result in equation (|38[) . the kernel -F lso (7) assumes the, still exact, form 

f iS °(7) = 4^^(7£(7)-X^T)-^ eH(7 .7 7q ' / ^ 

+2s 



M V wy ) " sin 2 a Jqr (<r 2 



, 9g(7 2 7qr) / dzu w (V3gK 2/s (g)Li /3 (g) - n) . (127) 
an a Jqr v ' 

The first term on the first line of equation (|127[) is the non-resonant isotropic kernel, calculated by integration over the full 
physical domain. The second term on the first line is an integral at given 7 over the quasi-resonant domain correcting for the 
fact that the non-resonant domain really does not extend over the full physical one. The second line is the quasi-resonant 
contribution proper. The variables a, x and g are functions of 7 and zu that may be found from appendix [A] and equation 
(|94[) . The QR integration over zu is between zu qr - and zu qr + defined in the caption of figure [6] 

As 7 approaches infinity, the interval [zu qr -, zu qr +] becomes more and more symmetrical about zero while g(+w,^f) 
and g{—w, 7) converge to each other. The second and third terms in equation (I127[) then approach the integral of an odd 
function over an interval symmetrical with respect to zu = 0, which then vanishes. As for the circular polarization absorption 
coefficient Kiv calculated in section 15721 higher order terms determine the QR contribution. Since only isotropic distributions 
are considered in this section, corrections caused by the anisotropy of the distribution function are absent. However, the 
interval [zu q , — , zu qr +] is slightly asymmetrical with respect to zu = 0, which gives rise to a first order "offset" correction. 
Moreover, a(+zu,^y) and a(— zu, 7) slightly differ when zu <C 7ttsin 2 a, as can be seen from equation (|Aip . and then g(+zu, 7) 
and g(—zu,j) differ by the quantity gi in equation (| 105[) . Thus the integrands in the QR integrals in equation (|127|) have 
a small even part that gives rise to a "parity" correction to these integrals. At large 7's, the offset Azuqi is given by the 
expansions in equations (|.12|l (fj3j, eventually giving for the offset correction ( appendix [K}: 

^off(7) = -^2. 7?t cosa + 8 S9 7«cosa (VSK 2/3 (1)L 1/3 (1) - tt) . (128) 

The parity corrections are evaluated by integrating the even part of the integrands in equation (|127|l over the symmetric interval 
[-3 1/3 (7ttsin 2 a) 2 / 3 , +3 1/ ^ 3 (7usin 2 a) 2 / 3 ]. The even part of the integrand on the second line of equation (|127[1 is obtained by 
Taylor-expanding in g the function in the parenthesis about go as in equation (|105|l , then performing the integration by using 
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g as the integration variable. To sufficient accuracy, g is related to vj by the lowest-order relation in equation (|J6[) . The 
minimum value g m of g, reached at vj = 0, is given by equation (|H4[) and approaches zero as 7 grows larger. The parity 
correction from the second line of equation (I127|l is thus evaluated taking g m = 0. The calculation of the parity correction 
from the first line is more straightforward. All calculations done (appendix [KJI, it is found that 



-irSq'yu cos a 



8, , , 4, 
- ln( 7 ) + - In 



21n(4) + -ln(3) 



4 + 



8^3 



dgg-^(y^gK 2ri (g)L iri (g)^ 



(129) 



An asymptotic expansion of F ls ° for large 7's is finally obtained by expanding the first term of equation (I127|l up to order 
7, integrating by parts the last term of equation (|129[1 . which partly simplifies with the second term of equation (| 128|) . and 
numerically calculating the coefficients of the terms of the expansion. This gives 

'4 . 



Flf (7) = 7rs 9 7 u cos a 



■In 



/jyv_\ 

V sin a J 



1.26072439 



(130) 



6.4 Isotropic kernel of the Faraday conversion coefficient in the low-frequency limit 

The integral in equation (|126p can be separated into a non-resonant and a quasi-resonant contribution. This is best achieved 
by writing the kernel J/ 1 " (7) as the sum of the contributions to it from equations (|119[) and (|120|l . which may be modified 
by extending the former integral to the full a-vj domain and adding to equation (|120fl a correction to account for the undue 
integration over the quasi-resonant domain so introduced. The variables are then changed to 7 and to. The extended non- 
resonant contribution then provides the result in equation (|39l) while the expression for the correction to the quasi-resonant 
contribution is similar to that in equation (|35f) . the boundaries of the integral over vo being however placed at ■uj qr+ ( , y) and 
ro gr _( 7 ) on the edge of the quasi-resonant domain shown in figure [5] The kernel H' so (7) then takes the form 

H (7) = - - sm a (7 (27 -3) Vf 3 ! + C(l))+ ~^ dw - [ j^^fi + -0 (g2 _ x2)7/2 ) 

Q H (j - 7 qr ) f^ir+ f 4x 2 /a - x\ 2 r . ,. r ,. 2u? 2 ( a - x \ „ ,. T 2vj 2 + a 2 . \ , lnl , 
+ \' 2 yq ' / dw -= \ K 2/3 ( g )L 2/3 (g) + —( )K 1/3 (g )L 1/3 (g - n °= . 131 

The first term on the first line of equation (|131[) is the non-resonant isotropic kernel, calculated by integration over the full 
physical domain. The second term on the first line is an integral over the quasi-resonant domain correcting for the fact that 
the non-resonant domain really does not extend over the full physical one. The second line is the quasi-resonant contribution 
proper. The low-frequency limit of this kernel is to be obtained from an asymptotic expansion in 7 of the right hand side of 
equation (| 13 1 f) . This low-frequency development is valid when 7 3> 7 qr ~ (it/(3 sin a)) 1 ^ 2 . The domain of validity of the low- 
frequency limit is where ui <<iC 3 7 2 sin a | fi| (equation (|27|) ~). which is equivalent to ui <§; aj c (7, a), where uj c is the characteristic 
frequency of synchrotron emission defined in equation (llOip . 

As discussed in section I5TT1 the variable g in equation Q94p is less than unity in the QR domain. At the order required by 
the calculations in this subsection, the lowest order approximation to the values of ■aj qr ±('y) is sufficient, namely 

tt qr ±(j) w ± (37 2 u 2 sin 4 a) 1/ ' 3 . (132) 
From equations (|A1I) - (|A2[) . the integrand on the first line of equation (|131|) can be written in terms of zu and 7, or G — 
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7« sin a/3, as: 



2z 4 2 i 2 (4a 2 + i 2 ) 162G 4 / 2ujcosa zu 2 sin 2 a + a\ x 2 



( a 2_ 3,2)5/2 ( CT 2_ 2-2)7/2 (^2 + ^2)5/2 ^ 3G QQ: 

405(JoG 4 / 2tucosQ to 2 sin 2 a + cr§ \ ( -. 2tt7Cosa 5vj 2 cos 2 a — U7 2 — a 2 \ , , 

+ (^+a 2 )7/2 + — 3G 9G 5 A 3G + 45G 5 J' ( 33) 

The primitive of this function is shown in Appendix [j] The integral over zu that appears in the second term of the first line 
of equation (|131[) can be obtained from it and expanded to the required order in 7. It is inferred from section \5A\ that this 
order should be 0(7 4/3 ). The expansion procedure, outlined in Appendix [j] leads to 

" qr+ f 2I 4 2 1 2 (4,J 2 + I 2 )\ .2 4 „ .2 2 7T . 2/3 / 7 U\ 

^t^t + o"n -r— s „, ^ ,» ~ 7r sin a 7 — 27r sin 07 sin ' a — 



dro ^ (^2. 3.2)5/2 +°o (g2 _ g 2)7/2 j « ^ sin a 7 - 27r sin a 7 --sin'a^— J . (134) 

The properly quasi-resonant contribution to H 130 ^) originates from the integration over zu of the second line in equation 
(|131[) . This contribution is calculated to order 7 4//3 by changing the variable zu for the variable g on which the K and L 
functions depend, much as was done in sect ion [5/fl but for the fact that the integration over zu is now performed at constant 7 
instead of at constant a. Some details are presented in Appendix[j] From these calculations it is found that the low frequency 
limit to the quasi-resonant contribution to the kernel _ff ls °(7) is 

HZnil) = 2 (7") 4/3 (sina) 2/3 3 1/6 £ dg g 2 ^ ^K 2/3 (g)L 2/3 (g) + K 1/3 (g)L 1/3 (g) - . (135) 



The integrand in equation (|135|l declines rapidly to zero out of the quasi-resonant domain (see section 15. 4[) and the value of 
the definite integral over g is, with a relative error of only 10 -4 , equal to 7r/(4 • 3 1//6 ). Adopting this value, we eventually get 

^ S r(7) « I (sina) 2/3 ( 7 ") 4/3 . (136) 

Accounting for the slightly different definitions of i? lso (7) and ^qr(c) in equations (|33p and ()121h . the results in equations 
(|136|) and (j 124f) are entirely equivalent since at fixed 7 the quasi-resonant domain is concentrated near zu — at a = 7tisin 2 a. 
This was not unexpected because the quasi-resonant domain is of a very small angular extent, so that the quasi-resonant 
contribution is insensitive to the fact that the distribution function is isotropic or otherwise. Finally, the asymptotic expansion 
in 7 of the first term in equation ()13ip . limited to the order 7 4//,i , is 



. 2 

— sin a 



( 7 (27 2 -3)v/7^T + /:(7)) ~-| sin 2 a(2 7 4 - 4 7 2 ). (137) 



The complete H ls °(-f) kernel in the low- frequency limit is the sum of the three contributions in equations (|134[) . (|136[) and 
(|137[) . The large terms proportional to 7 4 and 7 2 cancel out and the net result is, at low frequency: 



#lp(7) « £ (juf (sinaf [ 4 - -L ) . (138) 



g w > \ 1 y 34/3 

This result is in complete agreement with those in equation (|124l) and confirms that in the low frequency limit neither the 
resonant nor the non-resonant contribution to h predominates. 



6.5 Approximate expressions for the isotropic kernels over the full energy range 

Equations (|38l) - (|39p provide expressions for the isotropic kernels in the HF regimes that are valid in the interval [1, 7 qr ], 
although when 7 approaches 7 qr from below more terms should be retained in the Debye's expansion. Equations f) 130 fl and 
(|138p provide expressions of the kernels in the asymptotic LF regime at 7 ^ 7 qr . The transition between the HF regime and 
the asymptotic LF regime should take place over an interval 7 qr i < 7 < 7 qr 2, where 7 ql -i = Ai7 qr and 7 qr 2 = A27 qr and 
A\ < 1 and A2 > 1 are coefficients of order unity. For the kernel -F" sc (7) these coefficients should be close to unity because 
the change from the HF regime to the LF regime for this coefficient is dim. For H 1B °('y), or its derivative 7/' lsc (7), A 2 is 
more likely to be of order of several units since figure [7] shows that the asymptotic regime is in this case slowly reached. The 
figure [8] shows that exact kernels are continuous and at least once differentiable. Our proposed approximation adopts the 
HF expressions (|38|l and (|39[) at 7 < 7 qr i and extends them into the LF domain by a function that connects at 7 qr i to the 
HF expression and to its first order derivative and asymptotically merges into the LF limit in equations (| 130 p or (|138[1 . The 
kernels themselves or their derivatives could be interpolated that way, according to whether one is whishing to use equations 
(|33|l as they stand or to calculate them by integrating by parts. We opted for the second method for the coefficient h, which 
minimizes the inaccuracies due to the existence of singularities in the derivative of the distribution functions of truncated 
power laws. We interpolate from 7 ql -i to infinity by connection functions C c (7 — Tqri) that we define for coefficient c (= / or 
h) and element i in the connection (i = 1 for HF and i — 2 for LF) by 
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C ci (7 - 7qri) = I 1 - exp - 



(7-7qrl) 



(139) 

"ci /qrl / / 

The fact that this function is a square warrants that the interpolated function and its first derivative are continuous at 
7qri when the functions to be joined are. The extrapolations of the HF expressions to the LF domain that appear on the 
second lines of equations (|140p - (|141[) must then be continuous and derivable at 7 qr i but are otherwise unconstrained. A 
linear extrapolation proved satisfactory for F ISO , taking in this case 7 qr i = 7 qr 2 = 7 qr . Due to the non-negligible gap that 
the interpolation should bridge between the high-frequency and the low-frequency regime, the interpolation formula for the 
derivative _ff" so (7) of the kernel of the h coefficient is more sophisticated. A value of 7 qr i smaller than 7 qr has been chosen, 
and the extrapolation of the HF behaviour in the LF domain at 7 > 7 qr i has been endowed with a local bump that is meant 
to represent the cumulative effect of higher order terms in the Debye expansion when approaching its limit of validity. 

The decrement parameter \ C i is a constant, chosen to optimize the fit to the exact kernels in figure [8] For the kernel 
F lso we have adopted A/i = A/2 = 1/(1.7). For the derivative _ff" so (7) of the kernel of h we have adopted 7 qr i = 0.77 qr 
and different decrement parameters X C i in equations (|139[) , namely X^i = 0.4 and Am = 1-8. Other relevant parameters are 
described below. The resulting interpolation formulae for -F uso (7) and H' 1BO (-y) at 7 > 7 qr i are 

C f (l - 7qr) 7 (\ In ~ 1-26072439 ) + (1 - C,( 7 - 7ql )) x 

\ 3 V sm a / J 



*£?(7) = 



■nSnU cos a 



(47qr In ( 



7qr" 



Cw( 7 — 7qrl) 



7qr 



T / 2 ■ 

-— [u sma 
6 



-1-4 



2/3 



1% - 1 + 4 In ( 7qr + yffe - l) ( 7 - 7qr )) 



34/3 J 



1/3 



7r sm a 



4(27^-1)^-1 + 



4(6 7qr i - 5)7, 



qrl 



+ (1 - Chl{j - 7qrl)) X 

( 7-7qrl V 
\ M7qrl / 



( 7 - 7 qrl)+?7- 



1 + 



7~7qrl 
M7qrl 



(140) 



(141) 



The coefficient 77 = 103 and fi — A^i/1.5. The transition value 7qr rigorously is the value of 7 that causes the line of constant 7 
in figure[B]to tangent the NR-QR boundary. It slightly differs from the 7 value associated to a ql at vj = that is approximately 



representative of it and to which it reduces for cos a = 0. In our numerical evaluations we used the exact value, given by 



7 qr u sm 



(142) 



In appendix[L]we give some examples of Faraday coefficients derived from the interpolating formulae (|140[) - (|141[) and contrast 
them with the quasi exact ones obtained by double integration over zn and a of the expressions in equations (|25[) and (I26|) . 
The accuracy of the fit provided by the interpolating formulae is of order or better than about 10% in this frequency range. 



7 CONCLUSION 

7.1 Summary of the results in the different regimes 

We have derived an almost exact expression for the Faraday transfer coefficients in the form of equations (|25|) (|26l) . which we 
repeat here, as expressed in terms of the zu-a variables defined in equation (|23|) : 



/ 



h = 



n 2 2 o2 

2lV S„ 0J„ r \l 



m 3 c 3 , , dF 

5— azu da zu x — — 

sin a da 



■ dzuda 



da 



(x 2 j' a {x)K{x) - zu 2 J a {x)N a {x)) + 



OFo OFo 



dz 



do 



(143) 



(144) 



These expressions, which involve a two-real-variables integration, are completely general, applying to isotropic as well as to 
non-isotropic distributions and encompassing all regimes of particle-wave interaction. The Olver uniform expansions of high 
order Bessel functions presented in Appendix [G] can be used to speed up the integration in equations (|143[) and (|144p when 
very high orders a must be considered. In the high-frequency limit, when quasi-resonant contributions to these integrals are 
negligible, the transfer coefficients can be more simply written as 

_™2 



/hf 



u 2 n 2 

2tts -2— 



2 n2 



-dzu da 



dFp 

2 (a 2 -x 2 ) 3 / 2 J da ' 



/iHF — 



m 3 c 3 , , / 1 2x A {a 2 ~x 2 ) + alx 2 {Aa 2 +x 2 ) 



■ dzu da — 

8 



,2^/2 



dFp 
da ' 



(145) 
(146) 



Within the high-frequency limit, these expressions are also completely general. The integrands in equations (|145|l - (|146|l are 
more regular than those in equations (|143[) - (|146[1 . The integration over zu in equations (|145[) and (|146[) . although it is meant to 
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cover the non-resonant domain only, extends in the high-frequency limit to the full physical domain by lack of a quasi-resonant 
contribution, the QR domain then being ill-populated. 

We have particularized equations (|145[) - (|146|1 to a number of different physical situations. For isotropic distribution 
functions they reduce to equations (|38|l - (|39p . which coincide with known results for thermal distribution functions, as shown 
in section [3.2.41 A quadrupolar anisotropy, represented by a distribution functions as in equation (|5ip , produces the transfer 
coefficients shown in equations (|53|l - (|54|l while anisotropies of higher multipolar orders generate the coefficients compiled in 
Appendix |T] In all these cases, the expression of the transfer coefficients is reduced to a one-variable quadrature over the 
energy of the particles. For anisotropies that cannot be expanded in a sum of a few multipolar terms, the high-frequency 
coefficients are obtained by performing the double integrals in equations (|145|l - (|146p . either in these variables or in others, as 
shown in the case of a beam in section [3.3.21 

The high-frequency approximation progressively loses validity as the fraction of particles with a Lorentz factor 7 large 
enough to interact with the wave in the low-frequency mode increases. A regime change occurs when the quasi-resonant 
contribution ceases to be negligible, which, for a given frequency and a given direction of propagation, happens for Lorentz 
factors such that their characteristic synchrotron frequency, defined in equation (|10ip . is of order or larger than the frequency 
of the radiation considered. We offer clear physica l and mathematical explanat i ons fo r this behaviour, that has been previously 
well observed in the results of Shcherbakov ( 20081 ) and Huang fc Shcherbakov ( 2011 ) : the regime change occurs when particles 



in quasi-resonance make a contribution comparable to the non-resonant one. From a mathematical standpoint, the two- 
variable kernels in equations (|143p ~ (|144p must be represented differently in these two regimes, in terms of the Nicholson's 
approximations for quasi-resonance or in terms of the Debye expansion of high order Bessel functions for non-resonance. 

In full generality, the QR contributions to the transfer coefficients are given by equations (|99p - (|100p . The QR contribution 
to h is however better expressed as in equation (|120p . When quasi-resonant contributions are important, the NR and QR 
contributions should eventually be added and care should be taken to integrate the QR contributions in equations (|99p and 
(|120p over the QR domain only and the NR contributions in equations (|145fl — (|146|) over the NR domain only. We have 
compared the QR and NR contributions to the Faraday coefficients in section [5741 and shown that in the low-frequency limit 
none of them predominates. For a given frequency, the quasi-resonant contribution to the Faraday conversion coefficient h 
grows when a 3> o- qr to a value comparable to, and in fact numerically larger than, the non-resonant one. In the limit of 
large energies, the angular integration over the QR domain covers a very small interval in pitch angle space, as can be seen 
from equation (|9ip . Across this interval, the partial derivative d a Fo usually varies by only a little amount. Our estimation in 
section \5A\ of /iqr is for this reason expected to be a reliable one, that can be written in the low- frequency limit as 



^-da *-a^d a F ^ = Q,<T). (147) 

CUJ ' Gin /-v ) 



Equation (|147p is one among different contributions to h. It is is meant to apply only to values of a that are relevant to 
the asymptotic low- frequency limit. It is otherwise general and valid for any distribution function whatever the pitch angle 
distribution, the latter being anyway irrelevant to the quasi-resonant contribution. 

The complication of accounting for both NR and QR contributions has been overcome for isotropic distribution functions 
by actually calculating their sum in the low-frequency limit. This has lead to the results in equations (|33p . (|130p and (11381) . 
The asymptotic low-frequency contribution to the Faraday coefficients is in this case 

yft = ^f&Hi C os a r m 3 c 3 d 7 7 ( % In {^-\ - 1.26072439") ^° . (148) 
cuj 2 i 7 » 7qr \3 \smaJ J d-y 

2 l"l a/S ?/3ZT ^ 1 \ r „3s, ,/sdF 



cu** SW 8l 4 - 3^ J ^ ^ ■ ( 149 ) 

These simple analytical results, which apply to the asymp totic low-frequency regime 7 ^> 7q r , do not seem to have appeared 



inese simple analytical results, which apply to the asymp totic low-rrequency reg ime 7 ;g> 7q r , do not seem to have appearec 
so far in the literature. For example. IShcherbakovl (2008) and Illuang fc Shcherbakov! 1 20111 ) provide intermediate and low 



frequency regime fits to numerical results. For application to specific distributions functions, equations (|25p ~ (|26p , 1)31 [I (I32|) . 
(11451) (|146p or (|147p may of course be expressed in terms of any set of physical variables, such as the Lorentz factor 7 and 
pitch angle ■& of the particles, as was partly done in sections 13.21 and 13.31 

We proposed in equations (|140p - (|14ip approximate expressions of the kernels of the Faraday coefficients for isotropic 
distribution functions that interpolate between the high-frequency regime and the asymptotic low-frequency one. These 
approximations to the kernels generate reasonably accurate expressions of the coefficients / and h that can be written in the 
form of a simple one-variable quadrature over the entire energy range. In the high-frequency limit, we have similarly reduced 
to a one-variable quadrature on energy the expression of the Faraday coefficients for a large class of anisotropic distribution 
functions proportional to Legendre polynomials of low degree depending on the cosine of the pitch angle. These functions 
may be used as a basis for expanding distribution functions of any simple kind of anisotropy. For distribution functions with 
sharp anisotropies our results can be made completely explicit by performing a double integration over energy and angle. 
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7.2 Discussion and prospects 

Our analytical results have been deduced from a new, but straightforward, calculation of the anti-hermitian part of the 
conductivity tensor expressed in a form in which the usual sum over discrete synchrotron resonances is exactly replaced by an 
integration over a continuous variable. This particular representation of the conductivity tensor has not been hitherto used 
in this context. The integral over this continuous variable has principal value singularities that can be exactly reduced to the 
sum of a regular part and of a residual part that has been shown to be safely negligible. As a result, simple, general, quasi 
exact and non-singular expressions for the transfer coefficients have been obtained in a form that parallels that of the familiar 
expressions of the synchrotron emission and absorption coefficients, but differs from the form in which Faraday coefficients 
are usually expressed in the literature. Not-withstanding this difference our results exactly coincide with known exact ones 
for thermal distributions. 

Unlike the hermitian part of the conductivity, which describes synchrotron absorption, the anti-hermitian part is never 
dominated by the contribution of a small quasi-resonant population traveling close to the direction of propagation of the 
considered wave. On the contrary, non-resonant particles make most of the contribution to the non-dissipative coefficients in 
the high-frequency limit and cannot be ignored otherwise. 

The transfer properties of relativistic plasmas may sensibly differ from those of cold plasmas, which may be a useful 
diagnostic. For example the Faraday rotation coefficient / of a symmetric pair plasma, with equal densities of electrons and 
positrons having identical distribution functions, exactly vanishes but its Faraday conversion coefficient h does not. The study 
of a field-aligned beam in sect ion [3 . 3 . 2 1 also revealed that the Faraday rotation coefficient is enhanced in this case for radiation 
travelling in the beam. The knowledge of transport coefficients in all physical regimes, high frequency or low frequency, should 
allow optimal use of inversion algorithms of multiwavelength polarization data to yield the density and magnetic structure of 
the source and of the intervening medium. 

The reconstruction of the magnetic field distribution in volume, from multiwavelength pol arization observations of syn - 
chrotro n emitting astrophysic al media, has been attempted either partially ("Faraday synthesis", Brentiens fc de BruvnlEooi ) 



cnrotron emitting astropnysical media, nas been attemptea eitner partially (/Taraaay syntnesis , mentions iz ac cruvn ^uuoj 
or fully ( Thiebaut et ahlboioh . but in the latter case assuming a known distribution of the electronic population(s) responsible 



for the synchrotron emission and/or the Faraday rotation. In the case of cold plasmas (e.g. the interstellar medium of our 
Galaxy or nearby galaxies) where the circular polarization is negligible (both from emission and transfer), these assumptions 
about the underlying electr onic spatial distribut ion(s) are a real limitation to reconstructing the magnetic field structure. 
However, it was also shown ( Thiebaut et al. 2010h that in the case of relativistic sources where a single electronic distribution 



would be responsible for both synchrotron emission and radiative transfer effects (through their contribution to the dielectric 
properties of the plasma), the multiwavelength observation of circular polarization (in addition to intensity and linear polar- 
ization) should in principle allow the simultaneous reconstruction of the magnetic field and electronic spatial distributions. In 
this context, the Faraday rotation and conversion coefficients derived in this work are of particular interest. The diagnostic of 
relativistic jets from galactic nuclei and of pulsar winds by these methods would be of particular interest. It should be noted 
however that our present results only yield the transfer coefficients in the plasma rest frame and should be transformed to 
the observer's frame for use in the inversion algorithm. 
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APPENDIX A: VARIABLES SUITED TO THE RELATIVISTIC PARTICLE- WAVE INTERACTION 

The variables a and zu are defined in equation (|23|l . They depend on the propagation angle a of the radiation (figure [5} and 
may be substituted to other dynamical variables of a particle, such as p± and p« , or the modulus p of its momentum (or 
its Lorentz factor 7) and the pitch angle 1? of the particle's velocity. The variables a, zu are convenient since some often-met 
operators acting on the distribution function take a simple form when expressed in terms of them, as equation (|A5|I shows. 
Assuming the vacuum dispersion relation to be valid, equation (|23[1 relates zu and a to 7 and p\\/mc by 

P\\ \ ( P\\ \ CT — ZU COS OL P\\ (JCOSQ - ZU /.,n 

7 cos a — , a = u ( 7 — cos a , 7 = 5 , — — = 5 . (Al) 

mcl \ mc I u sin a mc u sin a 



From this, p± and the variable x defined in equation (|23[) may be deduced. Using the notation <7o = it sin a, defined in 
equation (|23l) . we have 
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2 2 2 2 

Pi cr — ZU — CTn , „ 

where <jo = it sm a . (A2J 

The physical region of the w-a plane is where a is positive and p\ positive, that is, a' 2 > zu 2 + a 2 . Thus a assumes a minimum 
value, ao that is usually large and can only be reached when zu — 0. The dynamical state of particles that are represented 
by points on the border of the physical domain are those for which p± = 0. For particles traveling at a velocity close to the 
speed of light, the sign of zu on the boundary of the domain is opposite to that of the field-aligned component of the particle's 
velocity, as can be seen from equations (|A1[) . Particles with vanishing zu are those for which 7 cos a = \J^j 2 — 1 cos 17, which 
means that the field-aligned particle velocity equals the field-aligned wave velocity. For large 7's this relation approximately 
reduces to § = a. When zu does not vanish, its sign and value are indicative of the direction of the velocity of the particle. 
This makes the variable zu akin to an angular one. By calculating the jacobian of the transformation from p±, pu to w, a, it 
is found that 

m 3 c 3 

p±dp±dpii = — 5 — 7 dzu da , (A3) 

v? sin a 

where 7 is the function of zu and a given in equation (|A1|) . The operators acting on the distribution function fo that are most 
often met in these calculations translate in zu-a variables as 

dfo u p±f dfo nna df \ dfo 1 f^dfo df \ 

"5 — = — — — ~5 h cos a — — I , — = • ( A4 ) 

op± m^c z 7 \ da dzu J dp\\ ■yrnc \ da dzu J 

n ,,v uv± f dfo df \ ( dfo nft\\ ojuv ± sin 2 a dfo . . 

D(fo) = — +cosQ— , [uj- — +fc||D(/ ) = — . (A5) 

mc V dzu da J \ dp± 1 mc A da 



Conversely, in terms of the Lorentz factor 7 and pitch angle ■&: 
dfo 1 / 9/o 1 / cos a 7 cos $ \ dfo 



da u sin 2 a \ dj sin ■& \ ^p~\ 7 2 — 1 / d-d 



(A6) 



dzu u sin a \ dj sin 1? \ ^72 _ 1 7 2 — 1 j dti 



dfo 1 / 9/o,l/ 1 7 cos i9 cos a \ dfo , 

COSQ i^- + ^n; ^ r 3 — 1 — hur ■ ( A7 ) 



APPENDIX B: THE CONTINUOUS-SPECTRUM REPRESENTATION OF THE CONDUCTIVITY 



This appendix gives a few hints to calculate the matrix elements G a b in equati on (I63C. th e func tions T a being defined in 
equation (|60[) . Similar calculations being described in some detail in the paper by Qin et al.l (|2007f ). we only sketch here the 
main methods that should be used. The definitions of some quantities in Qin et al. 1 20071 ') may however differ from ours by 
a sign and the gyrofrequencies are defined in their paper with a sign different from ours. Consider first the case in which 
a = b — 1, the matrix element then being 



Gn(a.x)= [ — [ exp ( i (aw + x(sm(<f> — w) — sin cj>)\ 

Jo 2yr Jo 2yr V v ' 



(Bl) 
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The identity sin(</> — w) — sm<j> = — 2sin(w/2) cos((/> — w/2) allows to integrate over <f) by using the identity 



d9_ 

2^ 



i y sin & 



= Mv). 



(B2) 



Jo is the Bessel function of index 0. For some other coefficients, the integration over <j) yields a result propo rtional to the first 
or second derivative of Jo. Changing the angle w into 2/3, Gii is reduced to the quadrature ( Watsonl (1922) chap. 5) 



Gu(cr,x) = l — e 2l<T/3 Jo (2a; sin /3) = e ltT7r J a (x)J- a {x) , 



(B3) 

Some other elements, like G31, involve after integrating over <f> the derivative of the Bessel functions Jo. The integration over 
w is then completed by integrating by parts. The calculation of elements G23, G33 and G22 is less straightforward because 
the integrals over w cannot be directly obtained from the identity in equation (|B3[) . This difficulty is overcome by reducing 
the matrix element to a form that can be integrated by parts, by using Bessel's equation for Jo(t): 

d 2 J (t) 1 dJ (t) 



+ 



+ J (t) = 0. 



dt 2 t dt 
Consider for example the matrix element 



dm 



dw 







2tt 



G 33 (<t,x) = j — cos</> / — e 

Changing <f> to 9 by c, 
^ dw 



cos 



w) e 



ix(sin(<l> — w) — sin < 



(B4) 



(B5) 



G 33 (a,x) 



9 + w/2 + n/2 and then w/2 to /3, G33 in equation (|B5|I is changed to 

<-2tt 



2tt 







d0_ 

2^ 



sin 



sin 



2 w 



2ix sin sin 



w d/3 



Jo (2x sin p) + sin 2 J (2x sin /3)) . (B6) 



On integrating over 9, use has been made of equation (|B2|) and of its second derivative with respect to y. The Bessel equation 
(|B4[) being written for an argument t — tp = 2x sin fi, the second derivative term with respect to t in equation (IB4|I is 
expressed in terms of derivatives with respect to fi. The first order derivative term proportional to dJo(tp)/d/3 happens to be 
tan 2 /? ( J' {tfi)/tf/). All terms proportional to J'o{tp) are then eliminated by using equation (|B4[) again, which yields 



4x 2 (jofe) + sin 2 p Jo(tp)) 



ri 2 Jofa) 
dp 2 



with 



£j3 = 2x sin p . 



Using this, the expression of G33 in equation (|B6|) becomes 



G33 = — 



dPe 2iafi d 2 J {tp) 



4:TYX 2 



dp 2 



dJo(2x sin P) 



Atyx 2 



dd 



+ 



2ia 
4x 2 



dP 2iaf3 dJo(2x sin/3) 
7T 6 dp 



(B7) 



(B8) 



which can be integrated by parts again, noting that the square bracket term in equation (|B8[) vanishes because Jo (t) vanishes 
at t = 0, to eventually give, using equation (|B3[) 



<*»M = [^^0 (2-in/?)] V J /"" f e 2 ^J {2xs^p) = ^ (e 2 ^ - l) + Je^fxJJ- 



d/3 



(B9) 



APPENDIX C: THE DISTRIBUTION ASSOCIATED WITH MULTIPLE RESONANCES 

Let us consider the limit for N — > 00 of the following function of a: 

f x _ e 2i,N,\ COSOTT 

V / sin (77r 

This function enters expressions involving the integration of its product with a regular function of a. When TV diverges, 
exp(2i7riV<T) oscillates infinitely rapidly, leaving in the limit a vanishing integral when multiplied by any regular function, 
which is the case in any interval containing no integer value. At integer values of cr, the denominator sina7r vanishes, as does 
also the factor (1 — exp(2iirNa)). To find the limit value, we expand by setting a = n + 77, where 77 is 0(1/N). This gives 

g(a) (1-e w g(n + 77) 2tt/V 77 - 1 2N — — — . (C2) 

V / sma-K \ 2-kNt] J 

Near the singularity, the real term, which is proportional to 2nN 2 ri, appears as a symmetrical cut replacing the diverging 
part of the function cos air/ sin o-k by an odd function of 77 that vanishes at 77 = in an interval that decreases with N as 
1/7V. This means that as N — > 00, the real part of the diverging integral of g(a)(cos air/ sin ait) should be understood as a 
principal value. The imaginary part is proportional near the resonance to 2N(sin2irNri/2-KNr]) which is large at 77 = 0, and 
then quickly decreases in the vicinity of a = n, oscillating rapidly on a scale 0(1/N). The integral of the product with a 
regular function g(a) in the vicinity of a = n can be evaluated by changing from 77 to y — 2nNrj, which yields, in this vicinity 
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"^)(- <!w! w) w - 4ff(B) i IZ s ^ir = - i9{n) - (C3) 



In the limit N — > oo, the function in equation (|C1|I then converges to the distribution 

lim ( 1 _ e ^N,\SE^^ V (S^) =P (S^l)- i7r y S ^a-n). (C4) 
iv-s-oo V / smcT7r \sm<77r/ \sin<T7r/ ^-^ 

n 

where So is a Dirac distribution. This leads to the result in equation (|69l) . The distribution D(cos an j sin o"7r) is causal, which 
was obvious from the beginning since the integration has been carried over positive delay times only. This explicitly shows up 
in the fact that when a is close to a singularity, at a — n say, the distribution T> reduces, according to Plemelj's formula, to 

l(v—-m5 D (<T-n)) =-- L— . (C5) 

The presence of the infinitesimal positive imaginary part at the denominator of the right hand side member reveals the causal 
character of the distribution (see the last paragraph of section 12.21) . 



APPENDIX D: THE TRANSVERSE COMPONENTS OF THE CONDUCTIVITY 



The components XYZ of the tensor M being given by equations (|68|) . its transverse components in the x y z frame are 
obtained from them by equations (|70[) - (|71[) . This first gives 

inV± ( dfo , „ A /„/COS<T7T\ w 2 . , ,i , „ , , a 

M x , x , = — i- f w _^° +fc || £>(/ ) ) (D(- 7 — — ) j' a (x)- J' a (x)K(x) + ^ 
M i i - SqnV± ' 

lvl x' y' 



*| \ V "P\\ x ) x \ "P-L yy \ \ sin (T7T / nxj 

My,x ' = { ~ s ^ ( sinQ rSh ~ COSQ n^) { u wi +Wfo) ) HSS) ^wJiw-^w^)-) . 

M y , y , = cos^^Ll^+^^/o))^^^)^^-^^^)-^ 
| S2* | \ op± J x* \ VsmCT7r/ an J 



o invn (( dfo „, , .(j^ /„/cos(T7r\ t2 , . T , , „ T , A k±D(fn) 

+ sin 2 Q — I w -^-fcx£> /o - MP J 2 (a: - J CT (x 7V CT (a: + W ° y 

S2* | \\ op 1 1 X J \ \sincr7r/ / 7ra 

fcifl(/o)-) (W^ZL) J*(z) - J a {x)N a {x) - — ) 

S2*| x \ dp\\ x J \ V sm an J an) 

inv\\ a ( df \ ( (cos<tk\ 2 1 

sm a cos a 



inv± a ( dfo 
sm a cos a — — — - — a> 



Although this is not obvious, the elements M x i y i and M y i x i differ by only a sign, since their factors of sin a are equal. This 
can be seen by calculating the difference of these factors of sin a, accounting for the definition of -D(/o) in equation (|20[) . The 
expression of M y i x i is improved by arranging the first parenthesis using the definitions in equation (|23[1 . The writing of M y i y i 
can be similarly improved by gathering the terms proportional to Bessel functions on the one hand, and those which are not 
on the other hand, then changing some factors accompanying sin 2 a and sin a cos a? by using equations (|20[) and (|23[) . Denoting 
D(fo) by D for brevity, this eventually gives, for example 

(»f£ + v)- <D1) 

- Ofr ^f) + 1 ("f: + ">") = 2 "" ("fr + k > D ) ■ < D2 > 

2 a f df \ . 2 fc ± D sin a cos a / f dfo a\ ( df 

—v± cos a — w- h fen X? ) +t)|i sm a 1 \v± \uj-^ k±D— +v\\ \ oj- h fcii-D 

x A \ dp± J x x \ \ dp 1 1 x J \ dpx 

._. , / lo sin a cos a dfo / , \ /Dsina cosa / dfo , /T ^„ N 

V fc ± oP|| V k±v± kj_v ± \ dpx J J J 

The last equation has been obtained by using the vacuum dispersion relation. The matrix elements can be given a simpler 
expression by using, instead of p± and py, the variables a and ru defined in equation (|23[) and presented in Appendix[X] This 
results in the equations (|72[) (|74l) . 
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APPENDIX E: RESIDUAL PRINCIPLE VALUE TERMS ARE NEGLIGIBLE 



The integrations which appear in equations 1)86[) and (|87[) are over the physical domain in the zu-a plane. This domain is 
bounded from below by the curve B of equation a 2 — zu 2 + u 2 sin 2 a. This constrains a to be larger than a minimum value, 
cro = usina, that is usually large (unless sin a is very small) because u = w/|f2| is large. For a given value of a, zu spans the 
interval [— zub{a), +zjt{a)], where zjt{a) = \J a" 2 — a' 2 . Each principal value term in equations (|86p and (|87[) can be written as 
a sum of a few contributions of the form 

' cos an \ . . . . dFo 



T = dzuda V 



^ ^ , x oro / ^,/cosa7r\ . ,,_,„. 
Oa l„ V sm an / 

(To 

where A(to, cr) is a function of a, w or x and B a (x) is the product of two functions that may be either J a (x) or its derivative 
J' a (x). The function C(a) is the integral over w of ^(ro, cr) B a (x) d a Fo on the interval [— h7{,(<t), +ro(,(cr)]. It varies on a scale 
of order a. On [n — 1/2, n + 1/2], C(cr) may be series-expanded by setting a = (a — n) and can be written, on this interval, as 

C(a) = C(n) + a C* (1) (n) + y C (2) (n) + • • • + ^ C* (fc) («) + ••• (E2) 



where C^ fc '(n) is the fc-th derivative of C(cr) at a = n. From this we get 



1/2 , /cosctttx , . ^2, t CW(n) , . f +1/2 , /a fc cosa7r\ , . 

do-PI- Ca=> A fc -i-i, where A fe = / da P — . E3 

n _i/2 Vsmo-Tr/ fc! 7_ 1/2 V sma7r / 

The integral defining Afe is regular when k > 0. The coefficient Afe vanishes when k is even because the integrand then is odd 
in a. Therefore the summation on k on the left of equation (|E3[) is on odd values k = 2m + 1 only. This series is absolutely 
convergent since each derivative C' 2m+1 ^(n) is expected to be of order C' 2m_1 ' (n)/n 2 , thus decreasing rapidly with m, so 
that the factorial (2m + 1)! at the denominator warrants a fast convergence. Moreover, h.2 m +\ also decreases with m. Let ni 
be the smallest integer for which m — 1/2 > cro, that is, the smallest integer for which the interval [ni — 1/2, ni + 1/2] is 
entirely contained in [cro, oo]. The expression T in equation (|E1 f) can be written as the sum of a contribution from the interval 
[cro j ni — 1 /2] and a sum over n > ni of terms similar to that in equation (IE3|l : 



/•ni-ir / x oo / oo r ,(2m + l) r \\ 

• /CT o n=nj \m=0 v ' / 

It is possible to replace the summation over n by an integral because the terms of the series in n do not change by much 
when n changes by unity. Noting that C ( - 2m+1 \a) is the derivative of C^ 2m '(cr), the second term can be exactly calculated by 
integrating over cr, which performs the summation over n as an integral, giving 

T= r\ ? M CM+ g A2m+1 ^ii ETfirst+Tsum . (E5 ) 

J \smairJ £-> n (2m + 1)! 



being the sum over all complete unit 
intervals. T sum is similar to the integral calculated in equation (|E3f) . for the particular value n = m, the derivatives c( 2m ) 
replacing the derivatives (7< 2m + 1 ). The ratio C^ 2m ' (ni)/C' 2m+1 ' (ni) is expected to be of order of m, which itself is close to 
cro. We may assume that there is an upper bound over all values of m, K cro say, to the ratio C' 2m ' (ni)/C' 2m+1 ^ (ni), so that 
we may evaluate an upper bound on T aum : 

C (2m+1) (m) _ „ r i+1/2 



T sum < Ka y A 2m+1 = Ka / da V\ C(a) . (E6) 

^ (2m +1)! y ni -i/ 2 \smairJ 



The contibution Tfl rst of the first interval, defined in equation ()E5I) . is itself of the same order of magnitude as the factor of 
K cro in the last term of equation (|E6[) . It then suffices to show that T sum , as estimated in equation (|B6[) . is negligibly small. 

Since ni is distant from cro by less than unity, the argument x of the Bessel functions that corresponds to a — ri\ < cro + 1 
is very small compared to its index ni ~ cto, the ratio x/a being of order ^/2/ao- Now, for arguments x y^cr, as is the case 



here ; the Bessel functions of the first kind J a {x) may be represented by their Carli ni approximation (|Abramowitz fc Stegun 



1964) that can be further simplified by using Stirling's formula for F(cr) at large a (jAbramowitz fc Stegunll 19641 ). This gives 



J a {x)^-= — , 3^)^-= — . (E7) 



/27rcr \ a J y/2ita 

The function A in equation ()E1|I may be either to 2 , zux or x 2 depending on the corresponding Bessel factor B a and the 
variable x is y/ a 2 — cr 2 — zu 2 . Setting zu = A \Ja 2 — a\, where A is in the interval [— 1, +1], we get 

/ A 2 (cr 2 -a 2 ) 



27TCT 



V 



A aja^-a 2 \ (l - A 2 )" ^ ( \^ a 2 - a 2 , a) , 
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where the upper element in the column vector corresponds to J 2 , the middle one to J a J' a an d the lower one to j' a 2 . The 
integral over A is itself a small number owing to the presence of the factor (1 — A 2 ) CT that is a large power of a number smaller 
than unity. But it is the factor in front of the integral that makes the whole expression extremely small, since a does not differ 
from oo by more than two units, implying that 



e 2 (a 2 



-I 



< 



4e 2 o- 



= O 



(E9) 



For frequencies at the peak of synchrotron emission, ao is of order of the square of the Lorentz factor 7 of the emitting 
particles, which means that the term on the right of equation (|E9|) is of order (7 2 ) ' , entirely negligible even at relatively 
small values of 7, such as 7 = 10, where its value is about 10 -200 . Then, setting a — ao + h where h is 0(1), substituting 00 to 
a wherever possible and otherwise approximating (a 2 — <tq) by 2aoh, the upper bound on T sum on the right of equation (IE6I I 
can be written as 

/ (2h/a ) 
(2h/a )? 

V 1 

Again, the prefactor proportional to (cro)~' TO causes the whole expression to be negligebly small. This results from the extreme 
smallness of the Bessel functions of the first kind and large indices when their argument is small. 



T < K 



■ns/2 



2e* 

(To 



(n-l-CTo) + i 



dh V 



(jH-ctq)- A 



h 



0-Q + -T 




\ 2 y^(\V2h^,v) . (ElO) 



APPENDIX F: NICHOLSON'S APPROXIMATION 



The Nicholson's approximation provid es expre s sions for the Bessel functions of large indices when their argument is close to 
the index, in our case slightly smaller. IWatson 1 19221 ). chap. 6, gives the following integral representations of these functions: 



Ja(x) — — cos (a9 — x sin 8) dd 



N^x) 



exp (— (at + x sinht)) dt . 



— / sin (a8 — x sin 9) dO H / exp (at — x sinh t) dt + 



exp (— (at + 2 sinh t)) dt . 



(Fl) 
(F2) 



For x and a large and almost equal, the integrand is rapidly oscillating, or decaying, and any of these integrals may be 
approximated by the method of stationary phase which consists in expanding the arguments of the trigonometric or exponential 
functions in the vicinity of the value of 9 or t where they are stationary. This gives, for example 



cos (a9 — x sin 9) d9 : 



((a - x)9 + x9 3 /6) d9 = 



7r .r 1 / 3 



cos 3 



2 1 / 3 (a 



32/3^1/3 



x , 3 

V + V 



dv . 



These integrals involve Airy-Hardy integrals of order 3 jPetiau 19551 ) defined by 



^3(9) = / cos(3gv + v ) dv . 



Sis(ff) 



sin(3(7t; + v ) dv , 



Ei 3 (<?) 



exp(— (3gv + v )) dv . 



(F3) 



(F4) 



The last, negative exponential, term in equations (|F1|I and (|F2|) being negligible compared to the others for large and positive 
x and a, we are left with 



J<j(x) 



31/3 

.TVS 



Ci 3 



2 l/3, 



32/3 3.1/3 



N a (x) 



71 X 1 / 3 



Si, 



2 1 / 3 (a 



32/3^1/3 



+ Ei 3 - 



2 1 ' /3 ( ( t 



32/3^1/3 



(F5) 



The function Ci3 reduces to a combination of modified Bessel functions of index ±1/3, but the functions Si3 and Ei3, 
separately, do not. However, the c ombination which appears in equation ()F5|) does. Correcting a little typo in Watson's book, 
p. 324 equation (8) (Watson 1922), we have 



n^/g cos I 
3 sin -| 



iff 



*•)) , sum + *(-,) - ■^ff> (/-j (**) + m (**)) 



(F6) 



The results in equation (|95|) are deduced from equation (|F6I) and from the definitions in equations (|93|l . Those in equation (1961 1 
result from the expressions in equation (|95|l by derivating with respect to x. Use is made of the recurrence relations for the 
functions I v , which translate into different recurrence relations for the functions K v and L„, that have different parity under 
a change of sign of the index 



K' k (y) = -K % (y)-±- K,(y), 



/„_i(y) + /„+i(l/) = 2I' v (y), 
L> k (y) = +L i{ y)-LL h (y). 



The results in equation (|96|) have been obtained by neglecting subdominant terms that appear as a result of the derivation. 
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APPENDIX G: OLVER UNIFORM ASYMPTOTIC EXPANSION OF BESSEL FUNCTIONS 



The Olver approximation for the Bessel functions and their derivatives at large order and for z < 1 is given bv lOlverl |l954) 



J a (az) 
N a {crz) 
j' a (az) 



1-z 2 



0-1/3 2^ CT 



+ 



r 5/3 



E 



4C \ 1/4 /Bi(a 2 / 3 C) ^ MO , Bi'(a 2 / 3 C) ^ B.(0 



1 - 



r l/3 



+ 



E 



CT 5/3 ^ ff * 
s=0 



4C \ - 1/4 / Ai(a 2 / 3 <) ^ C s (() , Ai'(a 2 / 3 C) ^ D S (C) 



z \ 1 — z 2 
4C 



r 4/3 



E 



7 2/3 



E 



Z \ 1 - 2 2 



-1/4 



Bi(a 2 / 3 C) ^ C.(0 , BiV /3 C)^^(0 



E 



7 4/3 fj" 
s=0 



,2/3 



(Gl) 
(G2) 
(G3) 
(G4) 



where Ai and Bi are Airy functions defined as in the appendix of the paper bv lOlverl (| 19541 ') and Ai' and Bi' their derivatives. 
The variable £ is related to z by 



3 2 

and A s , B s , C s and D„ are given by 



2s+l 



MO = E 6 ™r 3m/2 ^ 2s - m (^=) , c 1/2 s s (c) = -£a m <- 3m / 2 c/ 2s+1 _ m (- 7 =2==), 



m=0 
2s+l 



vT- «* 



C 1/2 C S (C) = -£& m C 3m/ V 2s+1 _ m (^=J==), d.(C) = E a m C" 3m/ V 2fl _ m (- ? =i 



m=0 V 1 Z ?7i=0 V 1 Z 

The Debye polynomials U n and V n are found by solving the following recurrence 

U n (x) = | (1 - a; 2 ) :rV n -i(x) + i jf (1 - 5t 2 ) I7 n -i(t) ctt, = x (x 2 - 1) (sECiC*) + ^U n -i{x)j + U n (x) , 

with Uo = 1 and Vb = 1 . The coefficients, 

(2s + l)(2s + 3)---(6s-l) L 6s + 1 
ao = l, 6o = l, a s = , b .=-— a, 

arise from asymptotic developments of the Airy functions at infinity, such as: 



(G5) 
(G6) 



A ^ /3 ^^7^-(-l^)E(-ir 



„-1/6a1/4 / o \ °° 

a s ..i, 2/3 an (T 1 (, l 11 3/2 \ v^, 

^W?- M ff '0~--^«p^--aC/Ji:(-l) 



CT 5£3s/2 



The Debye expansion is recovered by plugging these asymptotic expansions in equation (|G1|) - (|G4|) . e.g. 
1 



J a {az) 



■ exp 



2^3/2^ g I 



/2^(1 _ 22)1/4 

Making the following change of variables: 

z=-^T- e |c V2 =C-tanh^. 
cosh £ 3 

we recover equations (|109|1 and (|110p . In section T6.2I we make use of equations |l - (|G4|) to first order to compute numer- 
ically equation (I126|l . The corresponding implementation in Mathematica is available at the following URL: http://www. 
iap.fr/users/pichon/olver/. 



APPENDIX H: COMPARISON OF QUASI-RESONANT AND NON-RESONANT CONTRIBUTIONS 

We derive in this appendix the estimations presented in section 15.41 for the NR and QR contributions to the ro-integrated 
kernels of the Faraday rotation and conversion coefficients. Let us first consider the non-resonant contributions. Equation 
(|115[) gives /nr., which we may reduce to its dominant term, the first one in the parenthesis. Equation (|119[) gives /inr- The 
double integrals in these equations are over the non-resonant domain. The latter is represented in figure [6] and the equations 
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of its lower boundary B and upper boundary Sqr respectively are zu 2 + erg = a 2 and zu 2 + a 2 = 3 2//3 cr 4/ ' 3 . Taking advantage 
of the symmetry of the integrands, the ^-integrated kernels in equation (|121[) can be written as: 



9 f^°*° i-- 2 (*\ ^ 2 -(^ 2 + ^) mn 

K NK = 2 / (tow - , (HI) 



9 rf ^\ / 2(a 2 -(^ 2 + ag)) 2 ag(a 2 -( ro 2 + g g))(5a 2 -( S 7 2 + g g)) \ 

To evaluate these integrals in the limit of large er's, that is when a^/a — > 0, we substitute y = ro/cr to zu. The lower bound of 
the integrations becomes (3/cr) 1 ' 3 and the upper bound unity. The result is 



-(/) _ 

/(3/«r)V3 y "3 V 2 7 7 (3/CT) i/3 ~ a y 5 8V3 



We now turn to the quasi-resonant contributions which we evaluate in the same limit. Equation (|99|) gives /qr and equation 
()120p gives /iqr, combined with the boundary term /ibqr in equation (JTTSJ. The variable g is defined in equation (|94|) and 
can be expressed in terms of zu and a from equation (|A2|) . At a given cr, (/ reaches a minimum <? m (cr) at vj — 0. In the quasi- 
resonant domain, x may be considered equal to a except where the difference (a — x) is involved, which makes it possible to 
invert the relation between g and w, yielding, in the quasi-resonant domain: 

(ag + zu 2 ) 3/2 „2 , „2^/,„„2\ 2 /3 a-x _ 3 2/3 fg\V 3 ~ 3 



3 a 2 



2. 2 _ /„ 2\2/3 CT-K 3' fgy/" , . CT 

+ cr «(3fl>er) , rr^cr, « — - ( - 1 , g m (a )= — . (H4) 

X 2 Vcr/ 3cH 



Assuming as in section [5~4l that c^Fo depends linearly on zu, equation (|120|l can be recast in the form of equation (|12ip . When 
zu is changed for g and the limit a 3> o- qr ~ erf/ 2 is taken, the ^-integrated kernels assume the form 



K«l = V3v*f dgg^Jg^gJ 3 ( 



(«,)--£=) , (H5) 
K%1 = 2 (3 1/6 ) a 4 / 3 J 1 dg g 2 ' s (k, (g)L, (g) + Ki (jJL: (g) - -^j . (H6) 



The parenthesis in equation (|H5[) diverges as g approaches zero as <?~ 1//3 , but the integral still converges as g m approaches 
zero. Similarly, the parenthesis in equation (|H6[) diverges as g~ 4/ ' 3 at small values of g. Thanks to the factor g 2 ^ 3 , the integral 
over g is however convergent when its lower bound approaches zero, which it does when a grows very large. The integral 
in equation (|H5[) can be calculated numerically for <? m = 0, giving the result in equation (|123|) . From numerical calculation, 
it also appears that the right hand side of equation (|H6|I converges when a approaches infinity to a value that is close to 
(■7t/2) er 4//3 , with a relative accuracy of 10~ 4 . Adopting this convenient approximation, it is found that for <r 3> er 3 / 2 

<rH - \ - 4/3 , (H7) 

which is the result mentioned in equation (|124[) . For large values of g, that is out of the quasi-resonant domain, the parenthesis 
in equation (|H6[) decreases quickly to zero, owing to the fact that any product of a K and L function asymptotes to 7r/(v3<jf) 
in this limit. The remainder of the parenthesis decreases as g~ 3 whereas each of its terms considered individually decreases 
only as <7 _1 . The integrand thus declines with g as <? _7//3 , which would cause the integral in equation (|H6|) to converge if it 
were to be extended to a bound much larger than unity. Thus, the support of the integrand in equation (|H6|) does not extend 
out of the quasi-resonant domain. The decrease of the integrand in equation (|H5|) is however slower. 



APPENDIX I: ANISOTROPIC KERNELS OF HIGHER MULTIPOLAR ORDERS FOR FARADAY 
COEFFICIENTS AT HIGH-FREQUENCY 

Let us parametrize the distribution function as: F{pf, d) = Fn(7)Pn(cos After a bit of algebra, following the substitutions 
of Section 13.31 we find that the coefficient / can be written as in equation (|52[1 with the same C/ factor that appears in 
equation (|53[) . while the function D'j can be written as D'foil) + -D/iM^M- The factors D'f (left column) and D'f\ (right 
column) are listed below for n = 0, ••■ 6, where Pi denotes the value of a Legendre polynomials Pi at angle a, Pi = Pi (cos a): 
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4a/ 7 2 -1Pi 



7 

-6P 2 



4Pl 

2(2 7 2 +l)p 2 



7\/ 7 2 -l 

4(( 7 2 -l)p 1 -(ll 7 2 +4)p 3 ) 4((6 7 2 +9)p 3 +P 1 (l- 7 2 )) 

5 7 ^TT 5( 7 2 -l) 

54( 7 2 -l)p 2 -25(l0 7 2 +ll)P 4 6(-2 7 4 + 7 2 + l)p 2 +5(8 7 4 +24 7 2 +3)P 4 

21(7 2 -1) 77 ( 7 2_!)3/2 (II) 

8(ll7 4 -7 7 2 -4)p 3 -(274 7 4 +607 7 2 +64)p 5 8(2 7 4 + 7 2 -3)P 3 -5(8( 7 2 +5) 7 2 + 15)P 5 

18 7 ( 7 2 -l) 3/2 6( 7 2 -l) 2 

500(l0 7 4 + 7 2 -ll)P 4 -44l(28 7 4 +104 7 2 +33)P 6 20( 7 -l) 2 (8( 7 2 +3) 7 2 +3)P 4 -2l(l6 7 6 +120 7 4 +90 7 2 +5)P a 

660( 7 2 -l) 2 44 7 ( 7 2 -l) 5/2 

5( 7 2 -l)(274 7 4 +607 7 2 +64)p 5 -2(l452 7 6 +8132 7 4 +5175 7 2 +256)P 7 15( 7 2 -l)(8( 7 2 + 5) 7 2 +15)P 5 -14(l6 7 6 + 168 7 4 +210 7 2 +35)P 7 

130 7 ( 7 2 -l) 5/2 2 6 ( 7 2 -l) 3 



The first line corresponds to an isotropic distribution and reflects the result in equation (|38[) . The third line corresponds to 
the D'f part of equation (|53[) . Generically, we find that both coefficients scale like Q a [^y]P n -i(cosa) + Qp[y]Pn+i (cos a) with 
Q a and polynomials in 7. The dependence of the Faraday rotation coefficient on angle a for the considered anisotropy is 
set by these Legendre polynomials, P„-i and P n +i- The term independent of £(7) scales like l/(7 2 — l)'™ -1 '/ 2 and the term 
linear in £(7) scales like l/(7 2 — l) n//2 . Similarly, the function Df(-y) (equation l|52[)) can be written as Dfo(y) + D/i (7)^(7). 
The factors D/o (Jeft column) and D/i (right column) are gathered below for 71 = 1, • ■ • 6 (no Df term is present when n = 0). 
Again, Pi = Pi (cos a): 



2 f (57 2 +4)i 
37 \ 7 2 -l 



P 2 \ 2(-2 7 2 + (2 7 2 + 7)P 2 +2) 



-2 



3( 7 2-lp" 

(34 7 2 +86)p 3 -24( 7 2 -l)P! 6((2 7 4 + 15 7 2 +3)P 3 -2( 7 4 -l)p 1 ) 

5( 7 2 -l) a/2 5 7 ( 7 2 -l) 2 
12(-5 7 4 + 7 2 +4)p 2 + (74 7 4 +387 7 2 +64)p 4 3((8( 7 2 +13) 7 2 +63)p 4 -4(2 7 4 +5 7 2 -7)p 2 ) 

7 7 ( 7 2 -l) 2 7( 7 2 -l) 5/2 

(394 7 4 +35177 2 +1759)P 5 -20(l77 4 +26 7 2 -43)P 3 5((-87 6 -52 7 4 +48 7 2 + 12)p 3 + (8 7 6 + 160 7 4 + 195 7 2 +15)p 5 ) 

27( 7 2 -l) 5/2 9 7 ( 7 2 -l) 3 

6(276 7 6 +3764 7 4 +3789 7 2 +256)P 6 -20( 7 -l)( 7 +l)(74 7 4 +387 7 2 +64)p 4 15(2( -8 7 6 -96 7 4 +41 7 2 +63)P 4 + (l6 7 6 +456 7 4 +930 7 2 +215) P 6 ) 

88 7 ( 7 2 -l) 3 44( 7 2 -l) 7/2 

5(7 2 -l)(274 7 4 +607 7 2 +64)p 5 -2(l452 7 6 +8132 7 4 +5175 7 2 +256)P 7 15( 7 2 -1) (8( 7 2 +5) 7 2 + 15)p 5 -14(l6 7 6 +168 7 4 +210 7 2 +3s)p 7 
130 7 ( 7 2 -l) 5/2 26( 7 2 -l) 3 



(12) 



(13) 



The second line corresponds to the Df part of equation (|53|l . The same scalings are found w.r.t. P n _i(cosa) and P n +i(cosa). 
The terms independent of £(7) scale like l/(7 2 — l)' n+1 ' / ' 2 and the terms proportional to £(7) scale like l/(7 2 — i)("+ 2 )/ 2 . 
The anisotropic coefficient h can also be written as in equation (|52[) . with a prefactor Ch that is identical to that in equation 
(|54l) and a function D' h that can again be written as -D^oM + ^'h\ (7)^(7). The factors -D^ (^ e ft column) and (right 
column) are listed below for n = 0, ■ • ■ 3: 

iV7 2 -l(27 2 -3)(P 2 -l) 
(6 7 4 _17 7 2 -4)(P 1 -P 3 ) 
15 7 

7(4 7 4 -8 7 2 +19)+5(-20 7 4 +76 7 2 +3l)P 2 +36(2 7 4 -9 7 2 -8)P 4 
210V7 2 -1 

9((36( 7 2 -3) 7 2 +493) 7 2 +104)Pi+7(-132 7 6 +696 7 4 +509 7 2 -248)P 3 +100(6 7 6 -39 7 4 -80 7 2 +8)P 5 

1890 7 ( 7 2 -l) li^riyT 2 

The first line corresponds to the bracket in equation (|39[) and the third one to that in equation (|54[) . The terms independent of 
£(7) scale like l/(7 2 — 1)(" -1 '/ 2 while terms proportional to £(7) scale like l/(7 2 — l) n/ ^ 2 . The number of Legendre polynomials 
now increases with the anisotropy order n. Finally, the similar functions associated to Dh (defined as in equation (|52p 1 are, 
for n = 1, • • • 3: 

(2 7 2 + 13)(P!-P 3 ) 
15( 7 2 -l) 

7(-6 7 4 + 17 7 2 +4)+5(l8 7 4 +65 7 2 +4)p 2 -12(4 7 4 +37 7 2 +4)p 4 
105 7 ( 7 2 -l) 3/2 

9(64 7 4 -258 7 2 -33l)Pi-7(l68 7 4 +854 7 2 -197)P 3 + 100(6 7 4 +83 7 2 + 16)P 5 

630( 7 2 -l) 2 21U 7 ( 7 2 -l) 5/2 

The terms independent of £(7) scale like l/(7 2 — l)( n+1 '/ 2 and the terms proportional to £(7) scale like l/(7 2 — l)( n+2 '/ 2 . 



p 2 -i 
37 

P3-P1 

14 7 2 +94 7 2 P 2 -108 7 2 P 4 -7P 2 + 7 
42 7 -42 7 a 

100 7 2 P 5 -9(2 7 2 +3)Pi + (27-82 7 2 )P 3 



(14) 



(4 7 2 +l)(P 1 -P 3 ) 

5 7 ( 7 2 -l) 3/2 
(24 7 2 +5)P 2 -12(2 7 2 +l)p 4 +7 

500(4 7 2 +3) 7 2 P5+9(87 2 (7 2 -22)-7)Pi+7(-296 7 4 +127 2 +9)P 3 



(15) 
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APPENDIX J: FARADAY CONVERSION COEFFICIENT IN THE LOW-FREQUENCY LIMIT 

The equation g = 1 which defines the values zu qr ± of no at the edge of the QR domain for a given value of 7 can be exactly 
written as: 

/ \ 1/4 / \ 3/2 

W + <T 



f„2 , „.2\3/2 

( + er I /, so cos a 



37 2 M 2 sin a \ 7usin a 



1 + 



1 



V 



7 2 ^ 2 sin 4 a (l + rocos ^ Y J \ 2 2 \ 7 2 K 2 sin 4 Q 



1 1 



(Jl) 



7U sm- a 



Since | vj qr ± \ is much less than 711 sin 2 a, the three parentheses on the right are close to unity and the simplest approximation 
to vj qr ± in the low- frequency limit is given by equation (|132[) . Higher accuracy approximations to vj qr ± can be generated by 
successively iterating on equation f) J 1 1) . setting 7 -usin 2 a = 3G. The result, valid for | uj qr ± |2> Co, is: 

/„ 2 2 . 4 \i/3 /, 2cosa 1 / cos 2 a 5\ 1 / 10 cos 3 a 5 cosaA 1 \ ,,„, 

*v+ = +(37V sm * a ) / + + + m-)g+")> ( J2 ) 

/„ 2 2 . 4 \i/3 /"„ 2 cos a 1 / cos 2 a 5\ 1 / 10 cos 3 a 5 cos a \ 1 \ ,_„. 
= -(37-sma) _ ___ + ^_ _ _ j _ - (__ _ j _ + . . . j . (J 3 ) 

The primitive of the function in equation ()133[) is 



/ 2x 4 2 x 2 (4a 2 + x 2 )\ _ 81uj(7a 4 + Wzu 2 a% +4zu 4 ) ( G\ 4 36 cos a a 3 , (5o-g + 2ro 2 ) / G \ 3 

OT ( V (<7 2 - a: 2 )5/2 +<T ° ( a 2_ ^2)7/2 J - (C7 2 + a 2 )*/ 2 ( OT 2 + <7 2 )5/ 2 ) 

4sin 2 aro 4 + (9 - 7cos 2 a)cr 2 tu 2 + 5a 4 { G \ 2 2 sin 2 aow 4 + 5 sin 2 aaozu 2 + (3 - 2cos 2 a)oo / G 

-18zu ; — 5 , , — +12 cosa 



{vj 2 +a 2 f/ 2 V<W (^ 2 +<7 2 )5/ 2 \a 

(3 + 6 cos 2 a — 5 cos 4 o<)k7 4 +(6+14 sin 2 Qcos 2 a)<ToS7 2 +( 3 +8 cos 2 a — 6 cos 4 a)<r 4 



-2cos^(l + sin^a)ln^ ro 2 + a 2 + vjJ + vj^-^ '- v ' g + ^ )5/2 " ■ (J4) 

In the low-frequency limit UQ/\vj qr ±\ is small, as can be judged from equation (|132[) . this ratio being about (u/( 3 7 2 sin q)) 1//3 
which is meant to be small in this limit. The right hand side of equation (|J4[I can then be expanded in ao/tzr, and the definite 
integral over the interval [w qr -, zu qr +] calculated from it at any desired order in 7, using equations (|J2[) - (|J3[) . It will however 
be sufficient to extend the expansion to the first non- vanishing order in 7, which from section [5741 is expected to be 7 4 ^ 3 . The 
ratio (G/oo) is proportional to 7 and, as mentioned above, ao/zu qr ± is of order I/7 2 / 3 . The factors of all powers of (G/cto) 
in equation (|J4[) remain finite or approach zero as ao/zu approaches zero. This implies that terms proportionnal to (G/oo) 1 
and (G/cto)° are negligible at the order 7 4//3 , as are also the logarithmic terms. Similarly, the term proportionnal to (G/oo) 3 
does not contribute at this order because its factor is C(l/7 2 ). Those parts of equation (|J4|) that do contribute at order -f 4 ^ 3 
are the first term, the factor of which should be expanded to order (<jq/vj) 4 and the term proportional to (G/oo) 2 , the factor 
of which would have to be expanded to order (co/ro), had terms of that order been present in its expansion. This leads to 
equation (|J5[) below, which only involves the lowest order approximation to zu qr ± from equations (|J2[I - (|J3[) : 



x 2 (4cr 2 + x 2 )\ 81^7 / 1 o 4 \ ( G\ 4 72cc7. 2 /G x2 



\ 81 w 








I ~ M 
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VJ 4 J 





**> \ (-2 _ ,2 )5/ 2 + -0 {(j2 x2)7/2 — <**«{-) ■ (J5) 

We then obtain, using equation (|J5|) . the definition G = 7wsin 2 a/3 and the lowest order expressions for zu qr ±, the integral 
over w on the first line of equation (|131fl in the form of equation (|134l) . 



Let us now evaluate, along the lines described in section [U the quasi-resonant contribution to H ls °. In the low- frequency 
limit, the variable g assumes, for a given Lorentz factor 7, a value less than unity in the quasi-resonant domain, which reaches 
near zu — a, small minimum, g m , given by equation (|H4|) . The relation between the variables g and vj at a given 7 results 
from equations (|A1[) . (|A2|I and (|94[) . In the low- frequency limit, it becomes 

{■W 2 + (J 2 ) 3/2 2 , 2 2 - 4 \2/3 / 2/3 2/3\ / Tc x 

9~ , o ■ 4 . or ™ = (37 u sin a) [g ' - 9m } ■ ( J6 

3 7 2 it 2 sm a V / 

From equations (|J6|) and (|A2[) we also get 

) W (23/ 2 7 /sin 2 a) ' (J7) 
These results are used to convert the second line of equation (|131[l into functions of g, so that eventually: 

2tt 



^R(7) = -^(7«sin 2 a) 4/3 3 1/6 f dg (g 4/3 (K 2/3 (g)L 2/3 (g) + K 1/3 (g)L 1/3 (g) 



-^( 7lt sin 2 a ) 4 / 3 3 1 / 6 f d9 f<^V /S fjfi/3(s)ii/3(s) - ) ■ (J8) 
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The products K 2 / 3 (g)L 2 / 3 (g) and Ki/ 3 (g)L 1 / 3 (g) diverge at small g as g~ 4 ^ 3 - Therefore both integrals in equation (|J8[) 
converge in the limit g m — > 0, since the second integrand is less than g i ' 3 (Ki/ 3 L 1 / 3 — n/(\/3g)). Evaluating these integrals in 
the limit of vanishing g m , we obtain equation (|135[) of the main text. 



APPENDIX K: FARADAY ROTATION COEFFICIENT IN THE LOW-FREQUENCY LIMIT 

We derive here the low- frequency limit of equation (|127[) . which corresponds to 7 > 7 qr (equation (|98[) ~). The offset correction 
results from the small displacement Azu ql of the center of the interval [ro qr _ (7), , zu ql +(^)] with respect to zero. Approximations 
to ccqr±(7) are given in equations (|J2|) - (|J3|) . Azu ql being small, the integral of an odd function /odd (to) on [TO qr _, , z*7 qr +] 
is about twice Azu ql / O dd(+ro q ro), zu ql o being the dominant term in equation (|J2[). From the definition of the QR domain, 
g(-y,zu ql ±) is unity. This gives the offset correction associated with the the second and third terms in equation (|127[) . Since in 
the low- frequency limit (a — x) -C a and zu qr o <C o"o the raw expressions may be simplified to equation (| 128f) : 



Joff = ^ ^ n ^ q 7 ucosa + 8s 9 7mcosq (^/3K 2/3 (l)L 1/3 (l) - tyJ . (Kl) 

The parity correction is the integral over the symmetrical interval [— zu v o, +za ql o] of the almost odd integrands in equation 
(|127|) . Any of these two integrands can be factored as zu<&(zu). Only the small odd part of <&(zu), which is defined by 
4 > odd(7,^) = ($(7, +tu) — $(7, — CC7 ))/2, contributes to the integral over the interval [— T*7 qr o, +ra7 qr o]- Using the relation 
a 2 — x 2 — zu 2 + u 2 sin 2 a the function 5> dd associated with the integrand on the first line of equation (|127[) can be written as 

,m , . 1 4-yuzu sin 2 a cos a - . 

^(7, w) = - ' . (K2) 

A (zu 2 + u smr a) 

The associated integral in equation (|12Tf) can be calculated by integraton by parts, giving 



^° d zu - 2 - = in ( V^ ^^^ ! - 2 -=^= . (K3) 



'qrO ' ^qrO / y "VO 

This exact relation is expanded for large 7 up to terms of order 7 included, leading to an associated parity correction 

'8, , , 4 
3ln(7)+3 



^par = —irsq-yu cos a ( § In (7) + \ In ( ) + 2 In ( 4 .3 2/3 ) - 4 ) . (K4) 



go + gi where ffo(7, ro ) = ~ ~, . 2 \o and 9i = -2#o r-5- , (K5) 

1 rt 1 ofii ein 7t 



The parity correction associated with the integral on the second line of equation (|f 27[) involves the calculation of the odd 
part of the function gK 2 / 3 (g)L 1 / 3 (g), where 3(7, zu) is defined by equation (|94[) . For small zu/ju sin 2 a, g may be expanded 
by using equation (|A1|) and the second relation in equation ()A2|) into 

f (zu 2 + er, 2 ) 3/2 , to cos a 

q ~7 ■ 2 and f 1 = ~ 2 #° 

3 (■yusm a)' 1 711 sin a 

where gi <C go in the low-frequency limit. Denoting 3(7, ±tzr) by g±, the functions g±K2/3(g±)Li/s(g±) are then Taylor- 
expanded about go to provide an approximation of the odd part of gK 2 / 3 (g)L 1 / 3 (g), that is 

\ {g+K 2/3 (g+)L 1/3 (g + ) - g-K 2/3 (g-)L 1/3 (g-)) « _^£2|^ po _^_ ( g K 2/3 (g )L 1/3 (g ) j . (K6) 
2 71* sin a ago \ J 

The parity correction associated with the integral on the second line of equation 1|127[) can then be written as 
, {2) 8V3s,cosa f+™*ro d 
7usina J ago 



F (2) = _ gVo^a , ^ a flo ^ 2/3(5o)Li/3(5o) . (K7) 



Using equation (IK5|) . the integral in equation (|K7|I is expressed in terms of the integration variable go, noted g below, which 
ranges from g m = 170(7, w = 0) to unity. For large 7, g m — > 0. To the same order as in equation ()K1|) . it suffices to calculate 
the integral over g in equation (|K7|) to order 0(7°), which gives 



fiS = — 8\/3s„7ucos 



* f dg g-^ (gK 2/3 (g)L 1/3 (g)j . (K8) 

The corrections in equations (IK II) . (IK4[) and (|K8() are then added to the non-integrated term in equation (|127p . For consistency, 
the latter is expanded for large 7 as 7 In ^ \J^ 2 — 1 + 7 j — \/7 2 — 1 ~ 7 In (27) — 7. The integral over 3 in equation ()K8[) may 
be integrated by parts, resulting in the cancellation of the Bessel term in equation (|K1[) and leaving us with 



F\£(l) — KSq-yucosa 



1 in - 1 - 1 in(3) - 8 + ¥ r ^/aG^vaG?) * 



(K9) 



The integral over gi in equation (|K9|I is numerically evaluated to 2. f 62372 finally giving the result quoted in equation (|130[) 
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Figure LI. Left: quasi exact values of the Faraday rotation coefficient / for power-law distributions obtained by numerically calculating 
the double integral in equation l|25|l (full lines) compared to those obtained from the interpolating formula l|140|l for the kernel F lso (7) 
(dashed lines), for a range of values of the low-energy cutoff 7 m i n and for two different exponents n, as labeled. Right: Similar plots for 
the Faraday conversion coefficient h, from equations l|26ll and l|141[l . The radiation's parameters are u; = 15 | Q | and a = ir/4. 




Figure L2. Left; The thick line represents the ratio of the Faraday rotation coefficient f^ter f° r a thermal plasma as deduced from the 
interpolating formula for the ker nel in equation 11401 1 , to the high-frequency expression /gp deduced from equations J33D and 1381 . 
The thin line represents the fit bv lShc hcrbakov ( 200jJ) to the ratio of the exact value of / to t he high-frequency formula. These quantities 
are plotted as a function of the regime-change parameter X defined by [shchcrbakovl i2008l ) and in equation 14711 . Right: Similar plots 
for the thermal Faraday conversion coefficient h, from equations 11411 1 . 13311 and 139[l. The radiation' s parameters here are ui = 1000 | Q 
and a = 7r/4, allowing a comparison between these figures and figures 4 and 5 of IShcherbakovl ||2008T ). 



APPENDIX L: TRANSFER COEFFICIENTS FROM INTERPOLATED ISOTROPIC KERNELS 



The figures [Lll and [L2l illustrate how the interpolation formulae in equations (|140[) and 141 p perform in calculating the actual 
Faraday transfer coefficients in a plasma characterized by an isotropic power-law energy distribution with a low-energy cutoff 
(figure [Ll) or in a thermal plasma (figure |L2|) . Figure ILT1 compares the quasi exact values of the transfer coefficients in figure 
[3] to those derived from the interpolation formulae (1 14011 (I141|) for the kernel or its deri vative. The accurac y generally is of 
order 10% for h and somewhat better than 10% for /. Figure lL2l displays the ratio, which Shcherbakov 1 20081 ) refers to as the 
multiplier, of the Faraday transfer coefficients for a thermal plasma, calculated from the interpolated kernels (j 1401) and (114111 



Shcherbakov 



to the ir high-frequency expression deduced from equations (|38[) , (|39|l and 1)330 . In figures 4 and 5 of his paper _ 
( 20081 ) presents similar ratios of numerically computed exact transfer coefficients and corresponding high-frequency expressions, 
which he displays as a function of the regime-change parameter X defined by equation (|47)1 . He proposes fitting formulae for 



these ratios (his equations (28) and (29)). The figures 4 and 5 of IShcherbakov (|2008t ) can then be compared to our figures 
IL2I similar parameters having been adopted for the radiation. From this comparison it appears that within the range of 
parameters represented here our interpolation formula approaches the real value of the Faraday rotation coefficient somewhat 
closer than his proposed fit, while our interpolation also well represents the Faraday conversion coefficient, though slightly 
less accurately than his. 



